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SECTION  I 

COMPUTER  PROGRAM  DESCRIPTION 


INTRODUCTION 


The  programming  manual  describes  the  computer  program  and  analytical 
description  of  the  analysis  routines.  Section  I covers  all  aspects  of  the 
program  regarding  hardware,  software,  files,  system  implementation  and  the 
program  source. 


The  JOINT  program  is  essentially  a collection  of  composite  joint  analysis 
routines  that  have  been  interfaced  with  graphics  input/output  routines. 
Provisions  have  been  made  for  the  utilization  of  data  files  for  saving 
solutions  arid  printing  on  hardcopy. 
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OPERATIONAL  CONSIDERATIONS 

Hardware 

The  JOINT  program  utilizes  the  Tektronix  PL0T10  software  routines,  and 
may  be  executed  on  a hardware  that  can  emulate  the  Tektronix  4014/4015. 

The  screen  size  must  be  sufficient  to  handle  132  characters/line  and  64 
lines/page. 

The  4014  terminal  may  utilize  a joystick  or  thumbscrews  to  position  the 
screen  crosshairs.  If  there  is  a graphics  tablet  available,  the  user  is 
given  the  option  to  use  that  instead  of  the  crosshairs.  The  above  devices 
are  used  to  locate  and  transmit  screen  coordinates  for  the  selection  of 
various  options. 

Scope 

Joint  Types 

Bolted:  Balanced  Double-lap 

Unsupported  Single-lap 
Supported  Single-lap 
Stepped-! ap 
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Bended:  Double-lap 

Unsupported  Single-lap 
Supported  Single-lap 
Stepped-lap 
Scarfed 

Materials 

For  the  bolted  joints,  the  user  may  use  graphite-epoxy  materials  with 
either  25%  or  37.5%  zero-degree  plies.  Also,  the  user  may  choose  either 
steel  or  titanium  bolts. 

For  the  bonded  joints,  the  user  inputs  the  properties  of  the  adherend 
and  the  adhesive,  thereby  placing  no  restrictions  on  materials,  including 
metal . 

PROGRAM  CONSIDERATIONS 
Language 

The  source  was  developed  using  standard  FORTRAN. 

Standards 

CDC  standards  were  used.  The  availability  of  60-bit  words  made  the 
need  for  double-precision  unnecessary.  If  this  program  is  adapted  to  Install- 
ation with  a word  size  less  than  60  bits,  double-precision  accuracy  should 
be  used,  as  explained  in  the  analysis  routines.  > 

The  overlay  of  program  routines  was  accomplished  through  the  use  of  CDC 
OVERLAY.  Refer  to  figure  1 for  the  routine  overlay  map. 

Data  Files 

The  first  two  files  identified  on  the  CDC  program  card  are  for  the  SAVE 
and  PRINT  files,  respectively.  These  on-line  disk  files  have  default  local 
file  names  of  TAPE1  and  TAPE  2,  which  the  user  may  change  local  equivalencing. 
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OVERLAY  (JOINT, 0,0) 
PROGRAM  MAIN 
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Save  File 


The  SAVE  file  is  for  the  accumulation  of  user-designated  solutions  that 
may  be  used  as  Input  to  a problem,  or  written  t-  either  the  PRINT  file  or 
the  terminal  screen.  The  solution  Input  and  output  data  Is  contained  in 
the  program's  WORK  array.  When  the  user  electes  to  save  the  solution,  the 
program  positions  the  SAVE  file  at  the  end  of  data,  and  the  following  data 
Is  output  using  unformatted  write  statements. 

WRITE  (1)  NOT,  TMPNAM,  IT 

WRITE  (1)  NENT,  (WORK(I),  1=1,  NENT) 

where: 

i NOT  = sequential  solution  number  on  the  file 

TMPNAM  = 8-character  name 

| IT  = joint  type 

i NENT  = number  of  items  in  the  WORK  array 

WORK  = array  containing  solution  I/O  data 

The  end-of-data  on  the  save  file  is  indicated  by  NDT=999  in  the  first 
record. 

The  variable  ND  is  used  to  indicate  the  position  of  the  read/write  disk 
head.  ND  is  set  equal  to  the  solution  number  of  the  previous  I/O  operation. 

If  a read  of  the  save  file  is  indicated,  the  desired  solution  number  is 
compared  with  ND;  if  the  desired  number  is  not  greater  than  ND,  file  1 is 
rewound  before  the  search  is  made.  For  a save,  the  file  is  positioned  at 
the  end  of  data. 

The  ANAME  and  I TYPE  arrays  contain  the  name  and  type  of  all  solutions  on 
the  SAVE  file,  and  are  dimensioned  at  100.  Therefore  a 100  solution  limit 
to  the  number  of  solutions  contained  on  the  SAVE  file  is  due  to  these  array 
dimensions.  A consolidation  feature  is  available  to  the  user  to  purge  solu- 
tions from  the  SAVE  file.  All  except  the  indicated  solutions  to  be  purged 
are  copied  to  a temporary  file  then  back  to  unit  1.  The  solution  numbers  are 
resequenced,  and  corresponding  changes  are  made  to  the  ANAME  and  ITYPE  arrays. 
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The  user  may  elect  to  select  a solution  on  the  SAVE  file  to  use  as  input 
for  a problem.  The  ITYPE  array  is  searched  for  types  corresponding  to  the 
type  of  problem  selected,  and  only  those  corresponding  names  from  the  ANAME 
array  are  displayed  for  user  selection.  After  the  selection  is  made,  the 
SAVE  file  is  searched  for  that  solution,  and  the  data  read  into  the  WORK 
array . 

The  user  may  also  elect  to  output  selected  solutions  contained  on  the 
SAVE  file.  The  entire  ANAME  array  is  displayed  for  the  user  to  make  selec- 
tions. During  execution  the  SAVE  file  is  searched  for  a solution,  the  data 
is  read  to  the  WORK  array,  and  the  appropriate  output  routine  is  called; 
this  procedure  is  repeated  for  each  selected  solution  name. 

The  user  may  wish  to  input  a SAVE  file  that  was  constructed  from  a 
previous  execution  of  the  program.  A file  that  has  been  cataloged  on  CDC 
can  be  read  but  not  written  on.  Therefore  a copy  to  a local  file  is  usually 
required.  When  on  old  SAVE  file  is  input  to  the  program,  the  ANAME  and 
ITYPE  arrays  are  constructed  from  the  data  in  the  first  record  of  each  solu- 
tion on  the  file.  The  user  is  then  able  to  add  new  solutions,  read  input 
data  from  the  file,  and  print  from  the  file. 

Print  File 

The  PRINT  file  Is  for  the  formatted  output  of  solutions.  It  is  rewound 
and  created  new  each  session.  The  user  may  write  the  solution  to  a problem 
after  execution  and  display;  the  appropriate  output  routine  is  called,  and 
the  input  and  output  data  contained  in  the  WORK  array  is  written  to  the 
PRINT  file.  If  the  solution  is  output  to  the  SAVE  file,  the  user  may  use  i 

the  SELECTIVE  OUTPUT  mode  to  select  desired  solutions  to  be  output  to  the 
PRINT  file.  During  execution  each  solution  indicated  is  read  from  the  SAVE 
file  to  the  WORK  array  and  the  appropriate  output  routine  called  to  write 
the  solution  to  the  end  of  the  PRINT  file. 
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After  the  user's  session,  the  print  file  may  be  used  for  on-line  viewing, 
or  hardcopy  printing  via  a batch  job. 

Graphics 

In  addition  to  the  program  being  interactive,  it  is  executable  only  on 
the  Tektronix  4014/4015  hardware  noted  previously.  The  routines,  from 
PL0T10  Terminal  Control  System  User's  Manual  for  release  2. 0/3.0,  are  used 
to  plot  drawings,  change  character  size,  clear  the  page,  and  obtain  screen 
coordinates . 

Refer  to  the  System  Implementation  subsection  for  the  library  used  to 
supply  the  PL0T10  routines. 

The  use  of  both  Terminal  Control  System  (TCS)  commands  with  Fortran  I/O 
statements  must  be  handled  with  care.  TCS  commands  and  Fortran  WRITE's  use 
different  buffers  which  cannot  be  controlled.  Care  must  be  taken  to  follow 
these  guidelines:  ] 

o Before  calling  a TCS  command,  call  RECVR  if  any  WRITE's  precede 

it.  i 

1 

o Before  calling  Fortran  WRITE's,  call  ANMODE  if  TCS  commands  precede  j 

it.  | 

SYSTEM  IMPLEMENTATION  j 

1 

Update,  Compile,  and  Load  Procedure  j 

During  development  of  the  program,  the  procedure  outlined  in  Figure  2 1 

was  used  to  facilitate  updates.  The  three  separate  tasks  are  update,  compile,  j 

and  load.  The  absolute  file  JOINT  is  created  through  execution  of  this 

batch  job;  it  then  may  be  executed  on-line  (see  Figure  5).  \ 

» 

JOINTUPD  Is  the  update  file  used  for  the  CDC  UPDATE  command.  It  contains 
the  source  and  subsequent  updates.  The  system-produced  SOURCE  file  Is  not 
saved. 


J 
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JOINTLGO  is  the  LGO  file  output  from  COPYL.  It  contains  the  compilation 
of  all  previous  routines,  plus  those  compiled  In  the  current  run.  It  Is 
cataloged  only  because  It  Is  required  for  Input  to  COPYL. 

JOINT  Is  the  absolute  file  for  execution  on-line. 

When  new  routines  are  added  that  do  not  exist  on  the  JOINTLGO  file,  a 
full  compile  Is  executed  to  create  a new  JOINTLGO  file.  Figure  3 shows  a 
typical  setup  used  for  a full  source  compile.  This  batch  run  also  produces 
a new  JOINTUPD  file  from  the  source  file. 

Field  Length  Map 

Figure  4 shows  the  size  of  the  routines  In  the  program.  The  name  of 
each  subroutine  Is  Included. 

On-Line  Operation 

Figure  5 shows  the  procedure  to  execute  the  JOINT  program.  In  this 
example,  the  user  is  assigning  local  files  A and  B to  the  SAVE  and  PRINT 
files. 

SYSTEM  OVERVIEW 
Data  Flow  Description 

The  program  has  been  structured  so  the  WORK  array  Is  used  for  the  tempo- 
rary storage  of  all  Input  and  output  data  for  an  analysis.  The  Input  data 
Is  either  read  Into  the  WORK  array,  or  the  values  are  equlvalenced  to  WORK. 

The  WORK  array  Is  common  to  all  analysis  routines  where  the  Input  Is  received 
and  output  stored  In  the  WORK  array.  All  output  (OUT)  subroutines  use  data 
passed  In  the  WORK  array. 

If  re-analysis  Is  selected,  the  user  may  edit  entries  In  the  existing 
WORK  array,  or  copy  a SAVE  file  solution  to  the  WORK  array  for  edit. 

The  bonded  stepped-lap  analysis  routine  reuses  many  of  Its  variables. 
Therefore,  It  was  necessary  to  use  a temporary  WK  array  to  store  the  I/O  items. 


INSTRUCTION 


COMMENTS 
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ATTACH (OLOPL , JOINTUPO ) 

REQUEST (NEWFILE ,*PF) 
UPDATE(P,F,N=NEWFILE) 

RETURN (OLDPL .COMPILE) 
REWIND(NEWFILE) 

CATALOG(NEWFILE,JOINTUPD,RP=999) 
UPDATE (Q,P=NEWFILE) 

REWIND(COMPILE) 

REQUEST(LGO,*PF) 
FIN(I=COMPILE,OPT=2,R=3) 

RETURN (NEWFILE .COMPILE) 
ATTACH(OLDLGO.JOINTLGO) 

REQUEST (COPYLGO,*PF) 
REWIND(OLDLGO.COPYLGO) 

COPYL (OLDLGO  ,LGO  .COP YLGO ) 
CATAL0G(C0PYLG0,J0INTLG0,RP=999) 
RETURN (OLDL GO, LGO) 

REQUEST(ABS,*PF) 

REWIND (COP YLGO ,ABS) 

ATTACH(L,TEKLIB ,10=654321 ,SN=AFIT) 
LIBRARY(L) 

MAP (ON) 

LOAD(COPYLGO) 

NOGO(ABS) 

REWIND(ABS) 

1CATAL0G(ABS, JOINT, RP=999) 

*EOR 

*IDENT  MKS0109 


OLD  UPDATE  FILE 
FOR  NEW  UPDATE  FILE 
FULL  UPDATE  FOR  UPD  FILE 


COMPILE  ONLY  UPDATED  DECKS 

COMPILE  NEW  DECKS 
PREVIOUS  LGO  FILE 

REPLACE  OBSOLETE  DECKS 
NEW  LGO  FILE 

FOR  ABSOLUTE  JOINT  FILE 

LIBRARY  OF  TEKTRONIX  SUBR 

CREATE  ABSOLUTE  FILE 


****  DECK  OF  UPDATES  TO  FILE  **** 


*EOR 

♦COMPILE  0UT9, PRINT, ANAL! , BONDI ,ANAL9,ESCARF,PSCARF 


I Figure  2,  Update,  Compile,  and  Load  Procedure 
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INSTRUCTION 


COMMENTS 


ATTACH(OLDPL,JOINTUPD) 

UPDATE (P,F,S) 

RETURN ( COMPILE, OLDPL) 

REQUEST ( NEWFILE ,*PF) 

UPDATE ( F , I ^SOURCE ,N=NEWFILE) 

RETURN (SOURCE) 

REWIND(NEWFILE, COMPILE) 
CATAL0G(NEWFILE,J0INTUPD,RP=999) 
REQUEST (LGO,*PF) 
FIN(I=C0MPILE,0PT=2,B=LG0,R=3) 
RETURN( NEWFILE, COMPILE) 
CATAL0G(LG0,J0INTLG0,PP=999) 
REQUEST(ABS,*PF) 

REWIND(LGO.ABS) 

ATTACH (L ,TEKLI B , I D=654321 ,SN=AFIT ) 
LIBRARY (L) 

MAP(ON) 

LOAD(LGO) 

NOGO(ABS) 

REWIND (ABS) 

CATALOG(ABS , JOINT, RP=999) 

*EOR 

♦IDENT  MKS0203 


OLD  UPDATE  FILE 

FULL  UPDATE,  TO  SOURCE  FILE 

FOR  NEW  UPDATE  FILE 
OUTPUT  NEW  SOURCE  FILE 


FOR  NEW  LGO  FILE 
COMPILE  THE  INPUT  FILE 

NEW  LGO  FILE 

FOR  ABSOLUTE  JOINT  FILE 

LIBRARY  OF  TEKTRONIX  SUBR. 


CREATE  ABSOLUTE  FILE 


****  DECK  OF  UPDATES  TO  FILE  **** 

*EOF 


Figure  3.  Full  Compile  Procedure 


Figure  5.  Joint  Execution  Procedure 
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Subroutine  COPYWK  copies  the  WK  array  to  the  WORK  array  before  WK  Is 
reused  or  altered. 


I. 


i. 


Save  File 

All  data  In  the  WORK  array  Is  written  to  the  SAVE  file  In  subroutine  SAVE. 
The  number  of  entries  In  the  WORK  array,  NENT,  Is  set  In  the  analysis  Input 
routine.  Whenever  the  user  selects  to  SAVE  the  analysis,  the  file  Is  posi- 
tioned at  the  END  record  before  writing  the  data,  unformatted. 

Reading  data  from  the  SAVE  file  to  the  WORK  array  for  editing  Is  accom- 
plished In  subroutine  EDSEL  (EDI t SELection). 

Subroutine  PRINT  also  reads  the  SAVE  file  to  WORK  before  calling  the 
output  routines,  and  before  writing  to  TAPE  3 during  consolidation. 


Program  Routine  List 

The  following  is  an  alphabetical  list  of  the  program  routines. 


ANAL 

COPYWK 

0UT1 

READ1 

ANAL 11 

DBLB 

0UT4 

SAVE 

ANAL  7 

DESIGN 

OUT  7 

SELECT 

ANAL8 

EDSEL 

0UT8 

SID 

ANAL  9 

ESCARF 

0UT9 

STPLP 

B0LT1 

FKBOLT 

0UT11 

TH 

B0LT4 

FKPROP 

OPTION 

WDMAX 

BONDI 

FPROP 

PCT 

WDMIN 

B0ND2 

I NIT 

PCTSTP 

WT 

B0ND5 

MAIN 

PRINT 

WTSL 

BOX 

NAME 

PSCARF 

WWT 

BOXNO 

QUADMN 

XYLOC 

I 


I 
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Common  Blocks 

The  common  blocks  used  throughout  this  program  are  summarized  below  with 
the  variables  which  comprise  them. 


Unlabeled 

NOT,  NC,  NCT,  ITYPE(IOO) , IUNIT2 

/BLK1 / 

TMPNAM(2) , ANAME(100,2) 

/BLK2/ 

ND 

/BLK3/ 

NOCONV 

/BLOCK/ 

PR0P(6) , FBOLT (4) 

/CONSOL/ 

ICON 

/COPY/ 

KSTART(4) 

/DBLBLT/ 

MAT,  PR0P(6) , FBOLT (4 ) , TNOM,  DD,  MM,  FX,  ED,  FTUO, 
FSTN , LLF(9) , FFS(9) 

/I ANAL/ 

IA 

/TAB FLG/ 

I TAB 

/WK/ 

WK(800) 

/WORK/ 

W0RK( 2500) 

Variable  List  and  Descrii 


General 

Below  is  a list  and  description  of  frequently  used  variables.  Since 
the  variable  names  used  in  the  various  analysis  routines  are  unique, 
the  variables  and  their  equivalance  location  in  the  WORK  array  are  listed 
separately.  The  description  of  the  routines  covers  the  variables  particular 
to  each. 

ANAME  Array  of  analysis  names  contained  on  SAVE  file. 

IA  /I ANAL/  common  flag  for  Input  routines  which  are 

overlaid  by  analysis  routines. 

IC0N  /CONSOL/  common  flag  for  PRINT  routine. 
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IDVGT 


i 


ii 


i 


i 


I END 
IT 

I TAB 

ITYPE 

IUNIT2 

IX 

IY 

KC 

KD 

KT 

KN 

KSTART 

NC 

ND 

NOT 

NENT 

NU 

TMPNAM 

WK 

W0RK 


Flag  to  denote  divergent  analysis  for  bonded  stepped 
lap  solution. 

999  flag  signifying  end  of  data  on  SAVE  file. 
Analysis  type. 

/TABFLG/  flag  used  in  XYLOC  for  use  of  tablet  or 
crosshai r. 

Array  of  analysis  types  contained  on  SAVE  file. 

Set  = 1 to  print  message  upon  EXIT  that  PRINT 
file  contains  data. 

Integer  X screen  coordinate. 

Integer  Y screen  coordinate. 

Denotes  column  number  of  detected  name  on  screen 
(from  1 - 10). 

Denotes  row  number  of  detected  name  on  screen 
(from  1 - 10). 

Analysis  type  read  from  SAVE  file. 

Analysis  number. 

Array  of  starting  locations  within  WORK  for  bonded 
stepped-lap  data. 

Current  analysis  type  being  processed. 

/BLK2/last  analysis  number  processed  on  SAVE  file 
for  tracking  file  position. 

Total  number  of  solutions  on  SAVE  file. 

Number  of  entries  In  WORK  array. 

Unit  number. 

Temporary  contents  for  current  analysis  name. 

Temporary  work  array  for  bonded  stepped-lap 
routines. 

Array  for  all  Input  and  output  data  for  a given 
solution. 
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Analysis  Input/Output  Variables 


The  following  subsections  list  and  describe  the  various  I/O  variables 
used  for  each  analysis,  and  their  equivalent  location  In  the  WORK  array. 
The  number  in  parentheses  is  the  dimension. 


Bolted  Double-  and  Single-Lap  Joints 

variable 

work 

location 

description 

NX 

1 

input  joint  load  (lb. /in.) 

MATL 

2 

joint  material  code 

KB0LT 

3 

bolt  type  code 

FSTN 

4 

tension  failure  factor  of  safety 

TEMP 

5 

joint  temperature  (F.) 

W 

6 

bolt  spacing  within  a row  (in.) 

D 

7 

bolt  diameter  (in.) 

T 

8 

composite  plate  thickness 
(1/2  joint  thickness) 

M 

9 

number  of  bolt  rows  (output) 

WX 

10 

joint  weight  (lb. /in.) 

MM 

11 

absolute  value  of  M 

PB(9),  I PB (9) 

12 

decimal  and  % values  of  load  trans 

ferred 

FS(9) 

21 

margins  of  safety 

LF(9) 

30 

failure  mode  codes 

FX 

39 

computed  joint  load  if  NX  = 0 

MO 

40 

input  number  of  bolt  rows 

15 
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Bolted  Stepped-1 ap  Joint 


i. 


/ 


variable 

work 

location 

description 

M 

1 

number  of  steps 

FX 

2 

applied  joint  load 

TEMP 

3 

joint  temperature 

FSTN 

4 

tension  failure  factor  of  safety 

KB0LT 

5 

bolt  type  code 

MATH 

6 

inner  material  code 

MATL0 

7 

outer  material  code 

B(  9) 

8 

step  lengths 

D(9) 

17 

bolt  diameters 

W(9) 

26 

bolt  spacing 

WD(9) 

35 

W/D  ratios 

TI(9) 

44 

for 

each  i 

inner  material  thickness 

T0(9) 

53 

bolt  < 
row  | 

outer  material  thickness 

PB(9) 

62 

I 

% of  load  transferred  by  bolts 

PT 1(9) 

71 

% of  load  retained  by  inner  plate 

PT0(9) 

80 

% of  load  retained  by  outer  plate 

FS  (9) 

89 

margins  of  safety 

LF  ( 9 ) 

98 

failure  mode  codes 

F 

107 

critical  load  value 

IR0W 

108 

critical  row  number 

16 
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Bonded  Double-Lap  and  Supported  Single-Lap 


work 


variable 

location 

description 

K 

1 

load  type 

PL0AD 

2 

Input  applied  load 

0VRLAP 

3 

input  overlap  length 

(ADHESIVE  PROPERTIES) 

GAMMAX 

4 

max.  shear  strain 

ETA 

5 

thickness 

TEMP 

6 

operating  temperature 

TEMPO 

7 

cure  temperature 

TAUEL 

16 

elastic  shear  strength 

GEL 

17 

linear  elastic  shear  modulus 

TAUP 

18 

elastic-plastic  shear  strength 

GADHSV 

19 

non-linear  elastic  shear  modulud 

EPEEL 

8 

peel  modulus 

SGMAPL 

20 

peel  strength 

(ADHEREND  PROPERTIES) 

THICKI 

9 

thickness  (inner) 

THICK0 

21 

thickness  (outer) 

EINNER 

10 

Young's  modulus  (inner) 

E0UTER 

22 

Young's  modulus  (outer) 

GNUINR 

11 

Poisson's  ratio  (Inner) 

GNU0TR 

23 

Poisson's  ratio  (outer) 

ALPHAI 

12 

coeff.  of  thermal  expansion  (inner) 

ALPHA0 

24 

coeff.  of  thermal  expansion  (outer) 

FI 

13 

yield  strength  (Inner) 

F0 

25 

yield  strength  (outer) 

J 


variable 


work 

location 


description 


* 


l i 
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ETRNSI 

14 

transverse  modulus  (Inner) 

ETRNS0 

26 

transverse  modulus  (outer) 

SGTRNI 

15 

transverse  strength  (Inner) 

SGTRN0 

27 

transverse  strength  (outer) 

FCYINR 

28 

compressive  yield  strength  (Inner) 

FCY0TR 

31 

compressive  yield  strenghh  (outer) 

FSYINR 

29 

shear  yield  strength  (inner) 

FSY0TR 

32 

shear  yield  strength  (outer) 

FTYINR 

30 

tensile  yield  strength  (inner) 

FTY0TR 

33 

tensile  yield  strength  (outer) 

(OUTPUT  DATA) 
flfLAP 

34 

either  optimum  or  specified  overlap 

0VLAP 

35 

optimum  overlap  length 

GAMEL 

36 

linear  elastic  adhesive  shear  strain 

GAMMAE 

37 

non-linear  elastic  adhesive  shear 

strain 

GAMMAP 

38 

plastic  adhesive  shear  strain 

GAMMA 

39 

elastic-plastic  adhesive  shear  strain 

STRBND 

40 

elastic-plastic  adhesive  shear 
strength 

ELBND 

41 

linear  elastic  adhesive  shear  strength 

STRSAV 

42 

non-linear  elastic  adhesive  shear 

strength 

STRINR 

43 

Inner  adherend  strength 

STR0TR 

44 

outer  adherend  strength 

BNDPL 

45 

limit  load  due  to  adhesive  peel  or 

Interlaminar  adherend  tension 

18 


work 


variable 

location 

TAUMAX 

46 

IPRNT 

47 

JCRTND 

48 

ICRTND 

49 

FACT0R 

50 

TMPBND 

51 

Bonded  Unsupported  Single-Lap 

variable 

work 

location 

PL0AD 

1 

0LAP 

2 

(ADHESIVE) 

ETA 

3 

GAMMAX 

4 

TAUEL 

5 

GEL 

6 

TAUP 

7 

GADHSV 

8 

SGMAPL 

9 

EPEEL 

10 

(ADHEREND) 

THCKN 

11 

P0ISSN 

12 

SIGULT 

13 

description 

maximum  elastic  adhesive  shear 
stress 

stress  analysis  print  flag 

critical  end  for  strength  computation 

critical  end  for  stress  analysis 

number  of  bond  surfaces 

adhesive  shear  strength  for  given 
overlap 

description 

Input  applied  load 
Input  overlap  length 

thickness 

maximum  shear  strain 
elastic  shear  strength 
linear  elastic  shear  modulus 
elastic-plastic  shear  strength 
non-linear  elastic  shear  modulus 
peel  strength 
peel  modulus 

thickness 
Poisson's  ratio 
ultimate  yield  strength 
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variable 

work 

location 

description 

EY0UNG 

14 

Young's  modulus 

STRTRN 

15 

transverse  strength 

ETRNSV 

16 

transverse  modulus 

AKB 

17 

laminating  factor 

(OUTPUT  DATA) 

STR 

18 

remote  adherend  tension  strength 

BNDGLD(7) 

19 

combined  tension  + bending 

BNDSTE(7) 

26 

elastic  shear  strength  of  adhesive 

STRBND(7) 

33 

plastic  shear  strength  of  adhesive 

STRPL( 7) 

40 

limit  load  due  to  peel  or  Inter- 
laminar tension 

SIGAVG 

18 

ave.  applied  stress  remote  from 

joint  ■ 

•j 

SIGMAX( 7) 

19 

max.  induced  adherend  stress  at  edge 
of  overlap 

SIGBND(7) 

26 

peak  adhesive  shear  stress 

GAMBND(7) 

33 

peak  adhesive  shear  strain  j 

SIGBPL ( 7 ) 

40 

peak  adhesive  peel  stress  at  edge 

of  overlap  j 

1 

j 

It  Is  noted  that  the  last  two  output  groups  occupy  the  same  place  in  WORK. 

This  Is  so  because  only  one  group  of  variables  Is  computed,  depending  on  j 

PLOAD. 


Bonded  Scarf  Joint 


Since  the  same  routine  processes  the  Input,  refer  to  the  bonded  double- 
lap for  input  variable  descriptions.  The  only  differences  are  In  the  adherend 
variables  below: 


work 


variable 

location 

description 

T1 

9 

thickness  (left) 

T2 

21 

thickness  (right) 

El 

10 

Young's  modulus  (left) 

E2 

22 

Young's  modulus  (right) 

GNU1 , G1 

11 

Poisson's  ratio  (left) 

GNU2 , G2 

23 

Poisson's  ratio  (right) 

ALPHA1 

12 

coeff.  of  thermal  expansion  (left) 

ALPHA2 

24 

coeff.  of  thermal  expansion  (right) 

(OUTPUT  DATA) 

STRNG1 

34 

remote  strength  of  left  adherend 

STRNG2 

35 

remote  strength  of  right  adherend 

STR1 

36 

remote  stress  in  left  adherend 

STR2 

37 

remote  stress  in  right  adherend 

0VLAP(7) 

38 

array  of  overlap  lengths 

ELBND(7) 

45 

elastic  adhesive  shear  strength 

EEPBND(7) 

52 

transitional  adhesive  shear  strength 

EPBND(7) 

59 

plastic  adhesive  shear  strength 

BNDSTR(7) 

66 

peak  adhesive  shear  stress 

GAMBND(7) 

73 

peak  adhesive  shear  strain 

Kl(7) 

80 

critical  end  - elastic 

K2(7) 

87 

critical  end  - transitional 

K3(7) 

94 

critical  end  - plastic 

K4(7) 

101 

critical  end  - peak  stress 

SYM+2  108  code  designating  scarf  joint  for 

routine  B0ND1. 
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Bonded  Stepped-Lap  Variables 


* 


Due  to  the  variable  number  of  steps  possible  for  the  analyses  performed, 
and  the  re-use  of  the  various  arrays  after  each  analysis,  a temporary  WK 
array  is  used  to  equivalence  the  variables  during  analysis.  Below  is  a des- 
cription of  all  1/0  variables,  followed  by  tables  showing  their  equivalent 
location  in  WK,  and  relative  position  in  W0RK. 


KT 

no.  of  entries  used  in  W0RK 

SGNLD 

load  type  (-1 , 0,  1 ) 

NSTEPS .NCHECK 

number  of  steps 

DELTMP 

temperature  differential 

JTDBLR 

joint  or  doubler  flag 

NSYM 

flag  for  symmetrical  stress  distribution 

IFACTR 

single/double  bond  surface 

ETA 

adhesive  thickness 

TAUMAX 

peak  adhesive  shear  stress 

G 

elastic  adhesive  shear  modulus 

GAMMAX 

maximum  adhesive  shear  strain 

GAMMAE 

elastic  adhesive  shear  strain 

ALPHA0 

outer  adherend  coeff.  of  thermal 

expansion 

ALPHAI 

inner  adherend  coeff.  of  thermal 

expansion 

NOUT 

flag  for  C0PYWK 

IDVGT 

flag  denoting  divergent  analysis 

STEPL , STEP 

step  lengths 

TAU 

adhesive  shear  stresses 

GAMMA 

adhesive  shear  strains 

DELTA0 

outer  adherend  displacement 

DELTAI 

inner  adherend  displacement 

T0UTER 

outer  adherend  load 
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TINNER 


Inner  adherend  load 


THICK0,  THCKN0  outer  adherend  thickness 
THICKI,  THICKNJ  inner  adherend  thickness 


ET0TR 

ETINR 


outer  adherend  extensional  stiffness 
inner  adherend  extensional  stiffness 


STR0TR , STRGTR  outer  adherend  strength 
STRINR,  STRGNR  inner  adherend  strength 


TABLE  1 . BONDED  STEPPED-LAP  WK  ARRAY 


KT 

SGNLD 

NSTEPS 

DELTMP 

JTDBLR 

NSYM 

IFACTR 

ETA 

TAUMAX 

G 

GAMMAX 
GAMMA E 
ALPHA0 
ALPHAI 
N0UT 


STR0TR 

STRINR 


THICKI 

TAU 

GAMMA 


TINNER 

STRGTR 

STRGNR 


TABLE  2.  BONDED  STEPPED-LAP  WORK  ARRAY 


Base  values  not  dependent  on  step  data. 


INPUT  STEP  DATA 

7 arrays  dimensioned  at  NSTEPS  (a) 


ELASTIC  ANALYSIS 


11  arrays  dimensioned  at  NSTEPS  (b) 


ELASTIC-PLASTIC  ANALYSIS 


11  arrays  dimensioned  at  NCHECK  (c) 


TABLE  2.  BONDED  STEPPED-LAP  WORK  ARRAY  (Continued) 


1 

NCHECK 

(d) 

t 

IDVGT 

STEP 

1 

ELASTIC-PLASTIC  ANALYSIS  , 

THCKN0 

THCKNI 

TAU 

GAMMA 

DELTA0 

DELTAI 

T0UTER 

TINNER 

t 

WITH  INFINITE  ADHEREND  STRENGTH  ; 

9 arrays  dimensioned  at  NCHECK  (d)  J 

| 

1 

JOINT  PROGRAM  ROUTINE  DESCRIPTION 

The  following  is  a description  of  each  of  the  JOINT  routines  which  will 
include  the  following  items: 

1.  Algorithm  - an  overall  description  and  purpose  of  the  routine. 

2.  Argument  List  - a list  of  the  arguments  passed  to  or  from  the  calling 
routine. 

3.  Common  - a list  of  the  common  areas. 

4.  Length  - octal  words 

5.  Subroutines  Called  - a list  of  externals.  (PL0T10  routines  excluded) 

6.  Subroutines  Called  By  - a list  of  routines  that  call  this  routine. 

7.  Input/Output  - a description  of  reads  and  writes. 

8 Error  Handling  - a description  of  how  this  routine  handles  errors 

or  error  codes. 

9.  Flow  Chart  - a functional  diagram  of  the  main  portions  of  the  routine. 

10.  Symbol  List  - description  of  the  unique  variables  used  on  this  routine, 
not  covered  by  the  main  variable  list  in  the  System  Overview  subsection. 
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ANAL 


1.  Algorithm  - This  routine  analyzes  the  design  of  a bolted  double-lap 
joint. 

2.  Argument  List  - X,  F,  FU,  DFU,  NF 

3.  Common  - DBLBLT 

4.  Length  - 341  Octal 

5.  Subroutines  Called  - FKPROP,  FPROP,  PCT. 

6.  Subroutines  Called  by  - DBLB 

7.  Input/Output  - None 

8.  Error  Handling  - None 

9.  Flow  Chart  - None 

10.  Symbol  List  - 


A 

FSU 

GGMIN 

MATL 

BS 

FTU 

GMIN 

m 

D 

FTUO 

I 

NF 

DD 

FU 

J 

PB 

DFU 

FX 

JMIN 

PROP 

ED 

F2 

KT1 

T 

F 

F3 

KT2 

TNOM 

FBOLT 

F4 

LF 

W 

FBR 

G 

M 

WD 

FT 

FSTN 

GG 

MAT 

X 

ANAL11 

1.  Algorithm  - Bonded  Stepped-lap  analysis  routine  for  joints  and 
doublers. 

This  routine  analyzes  the  Input  in  the  WORK  and  WK  arrays  and  the 
output  data  passed  back  through  WK. 

There  are  two  main  analysis  sections.  The  elastic  analysis  is  from 
card  240  through  card  490;  the  elastic-plastic  analysis  is  from  card  550 
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to  card  960.  To  compute  the  potential  bond  strength  of  the  elastic- 
plastic  case,  it  Is  re-analyzed  with  the  Inner  and  outer  material  strengths 
set  very  high. 

The  routine  attempts  to  reverse  joint  geometry  automatically  when  the 
basic  load  carrying  assumptions  are  backwards. 

The  two  analysis  sections  use  the  same  variable  names  and  re-analysis 
with  reversed  ends  if  possible.  Therefore,  It  Is  necessary  to  use  the  WK 
work  array  to  store  data;  when  the  end  of  an  analysis  is  reached,  the 
routine  calls  COPYWK  to  oopythe  contents  of  the  WK  array  to  the  WORK  array 
for  final  disposition.  ANAL11  keeps  track  of  which  analysis  has  been  com- 
pleted with  N0UT. 

N0UT 

1 input  data 

2 elastic  solution 

3 elastic-plastic  solution 

4 El. -PI.  solution  for  bond  strength 

This  is  a self  contained  analysis  routine  that  has  been  adapted  from 
the  elastic  and  elastic-plastic  programs  contained  in  the  NASA  Technical 
Report  CR-1 12237. 

2.  Argument  List  - none 

3.  Common  - W0RK  and  WK 

4.  Length  - 2467  Octal 

5.  Subroutines  Called  - C0PYWK 

6.  Subroutines  Called  By  - DESIGN  (from  B0ND5) 

7.  I/O  - none 
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8.  Error  Handling  - When  the  time  interval  for  convergence  exceeds 
the  arbitrary  value  of  10  CP  seconds,  an  error  message, 

INPUT  ERROR.  A STEP  CONTAINS  AN  EXCESSIVELY  HIGH 
G VALUE,  AND/OR  AN  EXCESSIVELY  LOW  ET  VALUE. 

is  displayed,  the  routine  is  exited  with  N0UT  = 0. 

9.  Flow  Chart  - Figure  6. 

10.  Symbol  List  - Refer  to  the  I/O  variables  for  the  bonded  stepped- 
lap  routines.  All  others  are  unique  to  this  routine  and  are  contained 
in  a compilation  listing. 


t 

i 
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ANAL  7 
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1.  Algorithm  - Bonded  analysis  for  double-lap  and  supported  single- 
lap joints,  distinguished  by  the  number  of  bond  surfaces  factor,  ANB0ND. 

All  I/O  Is  passed  through  the  WORK  array.  The  re-use  of  the  var- 
iable STRSAV  has  necessitated  the  storage  of  the  original  value  in 
STRSA1 , W0RK(42) . 

This  routine  is  the  result  of  converting  the  non-dimensional Ized 
version  described  in  NASA  Technical  Report  CR-1 1 2235- 

Analysis  of  the  joint  automatically  computes  the  adhesive  and 
adherend  strengths  and  the  optimum  overlap.  If  an  overlap  Is  specified 
by  the  user,  the  bond  shear  strength  is  computed.  If  the  load  Is 
specified,  the  adhesive  shear  stress  and  strain  Is  computed  for  either 
the  specified  overlap,  or  the  optimum  if  the  zero  Is  specified. 

2.  Argument  List  - none. 

3.  Common  - WORK. 

4.  Length  - 1164  Octal. 

5.  Subroutines  Called  - none. 

6.  Subroutines  Called  By  - DESIGN  (from  BONDI). 

7.  1/0  - none. 

8.  Error  Handling  - If  the  joint  load  exceeds  the  bond  and  adherend 
strengths,  the  joint  is  overloaded,  and  IPRNT  - 4 is  a flag. 

9.  Flow  Diagram  - see  Figure  7. 

10.  Symbol  List  - The  1/0  symbols  are  covered  under  the  Routine  Variables 
for  bonded  double-lap  and  supported  single-lap  joints.  All  others  are 
unique  to  this  routine  and  contained  in  a compilation  listing. 
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Figure  7.  ANAL 7 flow  diagram 
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ANAL  8 


If  PLOAD  is  specified,  this  routine  calculates  the  associated  internal 
stresses;  otherwise,  the  loads  are  determined  at  which  various  allowables 
are  exceeded. 

If  an  overlap  is  specified  by  the  user,  calculations  are  based  on 
that  overlap;  otherwise,  calculations  are  based  on  overlaps  for  a range 
of  1/ t ratios  for  the  given  adherend  thickness,  t. 

2.  Argument  List  - none. 

3.  Common  - WORK. 

4.  Length  - 1216  Octal 


5.  Subroutines  Called  - none. 

6.  Subroutines  Called  By  - B0ND2. 

7.  1/0  - none. 

8.  Error  Handling  - Failure  to  converge  will  print  the  appropriate 
message  listing  the  iteration  loops  that  were  active.  The  flag  NOCONV 
is  set  to  1 and  returned  to  B0ND2. 

Asterisks  are  printed  signifying  strain  failure. 

9.  Flow  Oiagram  - see  Figure  8. 

10.  Symbol  List  - Reference  the  subsection  on  the  1/0  variables  for  the 


('  bonded  unsupported  single-lap  joint.  All  others  unique  to  this  routine 

l 

\ are  contained  in  the  listing  of  variables  for  a compilation. 


Initial 

parameters 


P 


Figure  8.  ANAL8  flow  diagram 
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ANAL  9 

1.  Algorithm  - This  routine  analyzes  a bonded  scarfed  joint  with  either 
two  (symmetrical)  or  one  (asymmetrical)  bond  lines. 

If  the  overlap  is  zero,  a range  of  7 i/t  ratios  are  used  for  calcu- 
lation. 

For  a specified  load,  maximum  adhesive  shear  stresses  and  strains 
are  calculated;  otherwise,  joint  strengths  are  calculated.  | 

If  the  extensional  stiffness  of  the  right  end  adherend  is  greater  \ 

than  the  left,  the  ends  are  reversed  and  the  user  is  so  advised.  ■ 

2.  Argument  List  - FACTOR 

3.  Common  - WORK 

4.  Length  - 712  Octal 

5.  Subroutines  Called  - ESCARF;  PSCARF 

6.  Subroutines  Called  By  - BONDI  | 

7.  1/0  - none. 

8.  Error  Handling  - If  a joint  is  overloaded,  output  values  are  set 
high  so  as  to  print  asterisks  indicating  failure  cases. 

9.  Flow  Diagram  - See  Figure  9. 

10.  Symbol  List  - Refer  to  the  Routine  Variables  for  the  bonded  scarf 
joint  for  1/0  variables.  The  remainder  are  unique  to  this  routine. 

A 

i 

I 
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Figure  9.  ANAL9  flow  diagram 
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B0LT1 


1.  Algorithm  - This  is  the  Input  interface  for  the  bolted  double-lap 
and  single-lap  joints.  All  input  Is  requested  from  the  uer  and  passed 
to  the  analysis  routines  in  the  WORK  array. 

The  SAVE  file  is  not  utilized  for  input  due  to  the  small  number  of 
constraints. 

2.  Argument  List  - none. 

3.  Common  - blank  (NOT,  NC);  BLK1  (TMPNAM) , BLK2  (ND),  WORK 

4.  Length  - 743  Octal. 

5.  Subroutines  Called  - DBLB,  NAME,  OPTION,  0UT1,  SELECT 

6.  Subroutines  Called  By  - DESIGN 

7.  I/O  - Displays  on  unit  6 the  header,  typical  drawing,  constraint 
descriptions,  and  prompting  messages.  Reads  from  unit  5 the  requested 
input  data  values. 

8.  Error  Handling  - none. 

9.  Flow  Diagram  - not  necessary. 

10.  Symbol  List  - The  following  is  a description  of  unique  variables  for 
this  routine.  Refer  to  the  previous  subsection  on  Routine  Variables  for 
others. 

DIS  overlap  length  of  2 inches  for  picture 
IDEL  screen  x width  of  option  boxes 
II  return  code  from  SELECT 
IX  starting  x coord,  of  first  box 

IY  starting  y coord,  of  first  box  (top) 

PT  (not  used  - argument  from  DBLB) 

TH1  .2  In.  outer  material  thickness  for  picture 

TH2  .5  in.  outer  material  thickness  for  picture 
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1.  Algorithm  - This  is  the  Input  interface  for  the  bolted  stepped-lap 
joint  option. 

An  existing  solution  on  the  SAVE  file  may  be  used  as  basic  Input 
for  the  WORK  array.  Editing  Is  accomplished  the  same  way  as  the 
RE-INPUT  option. 

This  routine  describes  the  Input  data  required,  then  prompts  the 
user  for  that  data  by  list-directed  reads.  Processing  will  not  continue 
until  the  read  list  has  been  satisfied.  After  all  input  Is  read  in,  the 
user  is  given  options  to  RETURN  to  DESIGN,  EXECUTE,  or  RE-INPUT  data. 

The  RE-INPUT  option  for  editing  prompts  the  user  for  viewing  and 
updating  the  problem.  All  Items  that  are  requested  must  be  entered, 
even  If  only  one  item  Is  different  from  the  original.  When  editing  is 
complete,  the  three  options  are  displayed  again. 

After  EXECUTE  has  been  selected,  the  STPLP  analysis  routine  is 
called.  Upon  completion  the  output  routine  is  called  and  selections 
processed  by  SELECT.  If  re-analysis  is  to  be  performed,  II  is  returned 
unchanged;  if  the  user  does  not  wish  to  re-analyze  a SAVE  file  solution 
or  the  existing  input  from  the  previous  analysis,  all  new  data  is 
requested. 

2.  Argument  List  - none. 

3.  Common  - blank,  BLK1 , WORK 

4.  Length  - 1132  Octal 

5.  Subroutines  Called  - EDSEL,  INIT,  NAME,  OPTION,  0UT4,  SELECT,  STPLP 

6.  Subroutines  Called  By  - DESIGN 

7.  1/0  - Displays  on  unit  6 the  header,  picture,  required  constraint 
description,  and  prompting  messages. 
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The  read  statements  are  all  list  directed  for  the  problem  con- 
straint data.  Answers  to  questions  must  be  the  integer  equivalent  of 
YES  and  NO  as  displayed.  This  simplifies  the  testing  of  valid  responses. 


8.  Error  Handling  - A maximum  of  nine  steps  is  allowed.  A message 
will  be  displayed  if  an  invalid  step  or  number  of  steps  is  entered,  and 
correction  requested. 

9.  Flow  Diagram  - see  Figure  10. 

10.  Symbol  List  - The  list  and  description  below  excludes  the  1/0 
variables  conmon  to  the  bolted  stepped-lap  routines,  described  in  the 
Routine  Variable  subsection. 


IDEL 

input  to  OPTION  as  box  width 
code  1 , 2 or  3 

; returned  as  option 

IE 

return  code  from  EDSEL 

II 

return  code  from  SELECT 

IREV 

set  - 1 to  review  input  data 

IX 

starting  X location  of  first 

option 

box 

IY 

starting  Y location  of  first 

option 

box 

MI 

temporary  storage  of  old  no. 

of  steps  (M) 
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BONDI 
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1.  Algorithm  - This  routine  is  the  input  interface  for  the  analysis  of 
bonded  double-lap,  supported  single-lap,  and  scarfed  joints. 

The  analyses  are  differentiated  by  setting  the  value  of  W0RK(50)  in 
DESIGN.  ANAL7  uses  W0RK(50)  for  the  number  of  bond  surfaces  (ANBOND),  a 
value  of  1 for  the  supported  single-lap,  and  a value  of  2 for  the  double- 
lap. 


For  the  scarf  joint,  DESIGN  sets  W0RK(50)  equal  to  3.  The  user  is 
then  asked  to  enter  a 1 or  2 for  the  number  of  bond  surfaces  for  the 
current  scarf  problem;  this  value,  designated  SYM,  is  passed  to  the  analy- 
sis routine,  ANAL9,  through  the  argument  list  to  flag  the  problem  as 
asymmetrical  (SYM  = 1.),  or  symmetrical  (SYM  =2.).  A constant  of  2 is 
added  to  SYM  to  give  W0RK(50)  a value  of  3 or  4 for  use  throughout  BONDI. 
To  simplify  the  analysis  and  output  variable  equivalencing,  W0RK(50)  is 
copied  to  W0RK(108)  just  prior  to  calling  ANAL9. 

If  the  variable  LOOP  is  set  equal  to  1,  this  signifies  that  data 
exists  in  the  WORK  array,  whether  from  a SAVE  file  solution  or  the 
previous  problem.  Therefore,  the  edit  option  is  made  available  to  the 
user. 


Except  for  the  scarf  joint,  a tension  load  type  requires  adhesive 
peel  and  adherend  transverse  properties.  When  all  constraints  are  input 
(option  = 1),  the  load  type  is  entered  first  and  is  used  to  determine  if 
prompting  is  required  for  the  above  properties.  If  option  2 is  used  to 
edit  existing  data,  the  current  values  are  displayed  before  allowing  edit- 
ing; therefore,  it  is  necessary  to  request  the  load  type  prior  to  dis- 
playing the  values,  so  those  properties  may  be  omitted  or  included  on 
the  screen. 

The  routine  displays  the  appropriate  heading,  analysis  name,  and 
picture  on  the  screen. 

If  WORK  already  contains  input  data  for  editing  (LOOP  = 1),  then  the 
routine  skips  forward  to  display  all  the  current  parameters  on  the  screen. 
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If  all  the  constraints,  or  parameters,  are  to  be  Input  (LOOP  = 0), 
then  the  routine  displays  and  requests  values  for  the  three  basic  joint 
constraints.  Then  all  the  required  adhesive  properties  are  displayed 
and  values  requested  for  each  by  moving  the  screen  cursor  under  the  VALUE 
heading  for  each  constraint.  When  complete,  the  process  is  repeated  for 
the  adherend  properties,  and  LOOP  Is  set  equal  to  1. 

At  this  point,  all  required  input  values  have  been  read  into  WORK. 
Boxes  are  drawn  beside  each  constraint  and  the  RETURN,  EXECUTE,  and 
RE-DISPLAY  options  are  displayed.  A screen  coordinate  is  requested  to 
modify  a constraint  or  process  one  of  the  above  options.  If  an  edit 
box  is  selected,  the  screen  cursor  is  moved  under  the  MOD.  column  for 
that  constraint,  allowing  the  key-in  of  a modified  value.  This  cycle 
of  requesting  a coordinate,  selecting  a small  box,  and  entering  a modi- 
fied value  is  continued  until  the  user  selects  one  of  the  bottom  option 
boxes. 

If  the  user  selects  RE -DISPLAY,  default  values  are  set,  the  screen 
Is  cleared  and  re-displayed  with  current  values.  The  boxes  are  dis- 
played and  screen  coordinate  requested. 

RETURN  branches  the  user  back  to  the  analysis  options  immediately. 

EXECUTE  sets  allowed  values  to  their  defaults,  clears  the  page  and, 
for  the  scarf  joint,  calls  ANAL9.  For  the  other  analyses,  overlay 
restrictions  required  returning  to  DESIGN  to  call  ANAL7;  upon  returning 
from  ANAL7,  BONDI  Is  re-entered  with  I A = 1 to  branch  to  where  it  left 
for  calling  the  output  routine  for  displaying  the  solution. 

If  re-analysis  Is  picked  In  SELECT,  the  page  is  cleared,  and  the 
program  branches  back  to  call  EDSEL. 

2.  Argument  List  - none 

3.  Common  - blank,  BLK1 , I ANAL , WORK 

4.  Length  - 2455  Octal 

5.  Subroutines  Called  - ANAL 9,  BOX,  B0XN0,  EDSEL,  NAME,  0UT7,  0UT9, 

SELECT,  XYLOC 
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6.  Subroutines  Called  By  - DESIGN 


7.  1/0  - Read  from  unit  5 all  required  Input  data;  write  to  unit  6 all 
input  data  values  and  problem  headings. 

8.  Error  Handling  - If  an  Invalid  option,  number  of  bond  surfaces,  or 
load  type  is  entered,  re-entry  is  requested. 

9.  Flow  Diagram  - see  Figure  11. 


10.  Symbol  List  - The  following  symbols  are  not  covered  by  the  general 
list  and  description. 

DIS  - length  of  picture  overlap 

I DEL  - delta  screen  width  for  square  boxes 

IDELY  - delta  Y screen  value  for  placing  cursor  IDELY  down 
from  JBOXTP  coordinate 

II  - return  code  from  SELECT 

IT  - return  code  from  EDSEL 

IX  - X screen  coord. 

I XX  - X screen  coord,  for  MOD.  cursor 

1X1  - X screen  coord,  for  left  side  VALUE  cursor 

1X2  - X screen  coord,  for  right  side  VALUE  cursor 

IY  - Y screen  coord. 

IYY  - Y screen  coord,  for  MOD.  cursor 

II  - starting  box  number  for  draws 

JBOXTP  - array  of  Y coord,  for  top  of  small  edit  boxes 
K - load  type 

KK  - temporary  values  for  new  load  type 

K1  - X value  for  left  side  of  left  box 

K2  - X value  for  right  side  of  left  box 

LOOP  - set  equal  to  1 of  WORK  is  full  of  input  data 

N,N2  - number  of  boxes  for  DO  parameter 

SYM  - number  of  bond  surfaces  for  scarf  joint 

TH1.TH2  - picture  display  thicknesses;  outer,  inner 


CALL  XYLOC 


PLACE 
DEFAULTS 
INTO  WORK 


T 


EXECUTE 
CALL  OUTPUT 
CALL  SELECT 


Figure  11.  BONDI  Flow  Diagram  (continued) 
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BOND? 

1.  Algorithm  - This  routine  Is  the  Interface  to  the  bonded  unsupported 
single-lap  joint  analysis. 

The  routine  Is  basically  similar  to  BONDI,  except  the  differences 
were  such  that  modifying  BONDI  would  not  be  as  efficient  as  creating  a 
new  overlay  program. 

Since  the  BONDI  and  B0ND2  routines  function  so  closely  the  same, 
this  routine  will  cover  only  the  differences. 

There  is  no  need  to  request  a type,  since  only  a tension  load  is 
analyzed,  and  the  number  of  parameters  required  is  constant. 

The  only  additional  property  required  is  the  adherend  laminating 
factor.  The  temperature  properties  for  the  adhesive  and  the  adherend 
thermal  coefficients  are  not  required  for  the  analysis.  Also,  the 
adherend  properties  provide  for  a balanced  joint. 

2.  Argument  List  - none 

3.  Common  - blank,  BIK1 , BLK3,  WORK 

4.  Length  - 1231  Octal 

Subroutines  Called  - ANAL 8,  BOX,  BOXNO,  EDSEL,  NAME,  0UT8,  SELECT, 

XYLOC 

Subroutines  Called  By  - DESIGN 
1/0  - same  as  BONDI 


5. 

6. 

7. 
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8.  Error  Handling  - If  an  invalid  option  is  keyed  in,  re-entry  is 
requested.  If  the  problem  will  not  converge,  ANAL8  displays  a message, 
and  returns  NOCONV  = 1.  The  program  branches  back  and  requests  an  option. 


i 


9.  Flow  Diagram  - see  Figure  12. 

10.  Symbol  List  - Refer  to  the  general  variable  list  and  BONDI  except 
for  the  following: 

IYLINI  - Y screen  coord,  of  top  line 

12  - No.  of  boxes  in  column 

K3  - X screen  coord,  of  left  side  of  right  box 

K4  - X screen  coord,  of  right  side  of  right  box 

NOCONV  /BLK3/  - flag  that  solution  will  not  converge 
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Figure  12.  B0ND2  Flow  Diagram 
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B0ND5 


1.  Algorithm  - This  is  the  input  interface  routine  for  the  bonded 
stepped-lap  analysis  of  joints  or  doublers.  All  input  data  values  are 
equivalenced  to  either  WORK  or  WK. 

This  is  a routine  that  inputs  the  constraints  by  groups;  basic, 
adhesive,  and  adherend.  If  editing  is  desired,  the  user  must  re-enter 
all  the  items  of  the  group  requested  for  update.  Prior  to  execution  the 
input  data  is  copied  from  the  temporary  WK  array  to  WORK  by  calling  COPYWK. 
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When  execution  is  selected,  the  analysis  routine  is  initiated  by 
returning  to  DESIGN  and  overlaying  B0ND5  with  ANAL11.  After  execution 
B0ND5  branches  to  display  the  analysis  summary. 


2.  Argument  List  - none 


3.  Common  - blank,  BLK1 , IANAL,  COPY,  WK,  WORK 

4.  Length  - 2254  Octal 

5.  Subroutines  Called  - COPYWK,  EDSEL,  INIT,  NAME,  OPTION,  0UT11, 

SELECT 

6.  Subroutines  Called  By  - DESIGN 

7.  1/0  - Writes  on  unit  6 the  display  messages  for  input  requests. 
Reads  from  unit  5 the  parameters. 

8.  Error  Handling  - An  error  message  will  be  printed  if  an  Invalid 
step  number  has  been  entered. 

9.  Flow  Diagram  - See  Figure  13. 
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10.  Symbol  List  - The  following  were  not  previously  described  under 
general  program  variables,  and  the  bonded  stepped- lap  WORK  array. 


I 

i 


I 

i 

I 


f' 

i 


A 

- actual  strength  value  for  stepped  lap  sunmary 

B 

- allowable  strength 

ID 

- IDVGT  flag 

I DEL 

- X 

IDUM 

- dummy  argument 

IE 

- EDSEL  return  code 

II 

- SELECT  return  code 

IP 

- start  location  in  WK  (computed) 

IS 

- Y screen  step  depth 

IX 

- X screen  coord. 

IY 

- Y screen  coord. 

11 

- set  = 1 for  RE-INPUT  to  skip  request  for  name  & 

JTDBLR 

KSTART/COPY/ 

- starting  locations  in  WORK 

N 

- No.  of  steps  in  WK  array  output  data 

NOUT 

- counter  for  KSTART  and  COPYWK 

NS 

- temporary  storage  for  old  NSTEPS  value 

R 

- ratio  of  A/B 

RIM 

- minimum  R value  for  inner  material 

R0M 

- minimum  R value  for  outer  material 

SI 

- minimum  inner  B value 

S0 

- minimum  outer  B value 

TI 

- minimum  Inner  A value 

TO 

- minimum  outer  A value 

BOX 

1.  Algorithm  - This  routine  draws  a square  box  on  the  screen  given  the 
top  left  corner  coordinates. 

2.  Argument  List  - IX,  IY,  IDEL 

i 

3.  Common  - none 

4.  Length  - 47  Octal 

j 5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - BONDI,  B0ND2 
j 7.  1/0  - draws  a box 

i 

8.  Error  Handling  - none 

9.  Flow  Diagram  - not  required 

10.  Symbol  List  - none 


BOXNO 


f 

t 


i 


i 
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1.  Algorithm  - This  routine  determines  the  box  number  detected,  if  any, 
for  editing.  The  routine  takes  the  input  Y coordinate,  N,  and  returns  the 
box  value.  If  N is  returned  as  zero,  the  IVAL  number  does  not  fall  within 
a box. 

2.  Argument  List  - NMAX,  ITOP,  I DEL,  N 

3.  Common  - none 

4.  Length  - 51  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - BONDI,  B0ND2 

7.  1/0  - none 

8.  Error  Handing  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List 

IDEL  - height  of  each  box 

ITP  - array  of  top  coordinates  of  the  boxes 

IVAL  - value  of  N input 

IYDIFF  - difference  between  middle  box  top  and  IVAL 

N - Input  = Y coordinate 

output  = box  number 

NB  - bottom  box  number 

NMAX  - number  of  boxes 

NMID  - middle  box  number 

NT  - top  box  number 
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1.  Algorithm  - This  routine  copies  items  from  the  WK  array  into  WORK.  If 
J is  negative,  items  are  copied  from  WORK  to  WK  for  editing  in  B0ND5. 


The  KSTART  array  keeps  track  of  the  starting  location  of  the  data 
in  WORK. 

The  WK  array  contains  the  data  for  several  arrays;  N ITEMS  indicates 
the  number  of  arrays  in  WK  for  the  input  value  of  J.  IWS  is  simply  the 
starting  location  within  WK  for  each  if  the  items. 

The  simple  method  would  be  to  copy  all  items  of  WK  to  WORK,  but 
depending  on  the  number  of  items  in  each  array,  MAX,  much  space  could  be 
wasted  on  the  SAVE  file  when  WORK  is  copied. 

W0RK(1)  contains  a running  count  of  the  number  of  items  in  the  WORK 


array. 

2. 

Argument  List  - J,  MAX,  IDVGT, 

WK 

3. 

Common  - COPY,  WORK 

4. 

Length  - 145  Octal 

5. 

Subroutines  Called  - none 

6. 

Subroutines  Called  By  - B0ND5, 

AN ALII 

7. 

1/0  - none 

8. 

Error  Handling  - none 

9. 

Flow  Diagram  - none 
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Symbol  List: 


I - absolute  value  of  J 

IDVGT  - divergent  solution  flag 

IL  - last  item  in  array 

IS  - start  item  in  array 

IWS  - array  of  starting  locations  for  WK 

J - flag  denoting  type  of  data 

K - loop  counter 

KSTART  - array  of  starting  locations  within  WORK 
KT  - counter  for  number  of  items  in  WORK 
L - loop  counter 

MAX  - no.  of  steps 

NI  - no.  of  items 

NITEMS  - array  of  number  of  items  in  IWS  array 


- 


DBLB 

1.  Algorithm  - Search  routine  for  bolted  double-lap  joint  to  determine 
; minimum  wt.  i 

2.  Argument  List  - NX,  FST,  TCMP,  MATL,  KBOLT,  W,  D,  T,  WX,  M,  PB,  PT,  ' 

, LF,  FS 

3.  Common  - LBLT 

4.  Length  - 573  Octal 

i 

5.  Subroutines  Called  - ANAL,  FKBOLT , FPROP,  PCT,  QUADMN , TH,  WDMAX , 

WDMIN,  WT,  WTSL , WWT 

j 6.  Subroutines  Called  By  - ANAL! 

i 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Chart  - none 

10.  Symbol  List  - none 

.5 

J 
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DESIGN 


r! 


>■ 
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1.  Algorithm  - This  routine  displays  the  types  of  joints  available, 
prompts  the  user  for  a selection  and  branches  to  the  appropriate  routine. 

2.  Argument  List  - none 

3.  Common  - blank,  I ANAL,  WORK  I 

J 

j 

4.  Length  - 333  Octal  j 

5.  Subroutines  Called  - OVERLAY  j 

1 

6.  Subroutines  Called  By  - MAIN  j 

7.  I/O  - Displays  analysis  option  descriptions  and  reads  option  number.  \ 

i 

I 

8.  Error  Handling  - Checks  for  a valid  option  number  i 

9.  Flow  Diagram  - none  j 

10.  Symbol  List  j 

NGO  - option  number  j 
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I EDSEL 

n 

[I  1.  Algorithm  - This  routine  allows  the  user  to  select  a solution  from 

the  SAVE  file  for  editing. 

If  there  are  no  solutions  on  the  SAVE  file  with  a type  the  same  as 
the  current  analysis  type,  a message  stating  such  is  displayed.  Otherwise, 
all  the  names  with  a corresponding  type  are  displayed  on  the  screen. 

If  the  user  selects  RETURN,  IT  is  set  equal  to  zero  indicating  no 
| solution  from  the  SAVE  file  was  read. 

Otherwise,  the  user  selects  a name  by  its  screen  coordinate.  When 
a valid  name  has  been  picked,  the  screen  coordinate  is  converted  to  a 
! design  number,  KN,  and  the  solution  is  read  into  WORK. 

i 

2.  Argument  List  (IT) 

3.  Common  - blank,  BLK1 , BLK2,  WORK 

4.  Length  - 1001  Octal 

5.  Subroutines  Called  - INIT,  READ1 , XYLOC 

6.  Subroutines  Called  By  - B0LT4,  BONDI,  B0ND2,  BONDS 


[ 7.  1/0  - The  names  of  solutions  are  displayed  on  the  screen. 

[ After  the  SAVE  file  has  been  positioned  by  READ1,  the  data  is  read 


* 

f.  from  file  1. 


I 


10.  Symbol  List  - The  following  variables  are  not  Included  In  the 
general  list: 

IT  - set  equal  to  1 If  a name  picked 

IXMIN  - min.  X value  of  a name  on  screen 

IXP  - array  of  starting  X locations  for  up  to  10  names  per  row 

IYMAX  - max.  Y value  of  a name  on  screen 

IYMIN  - min.  Y value  of  a name  on  screen  computed  by  10  names 

per  row 

IYP  - array  of  Y screen  coordinates  for  each  row  of  names 

KC  - column  no.  of  name  (from  1-10) 

KD  - row  number  of  name  (from  1-10) 

KN  - design  number  of  solution 

TNAME  - copy  of  ANAME  with  blanks  for  names  with  non-correspond- 

ing types 
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Figure  14.  EDSEL  Flow  Diagram 
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ESCARF 


1.  Algorithm  - Elastic  Analysis  of  shear  In  unbalanced  scarf  joints. 

2.  Argument  List  - OL,  ETR,  THERMA,  LCRTND,  TAUAV6 

3.  Common  - none 

4.  Length  - 370  Octal 

5.  Subroutine  Called  - none 

6.  Subroutines  Called  By  - ANAL9 

7.  1/0  - none 

8.  Error  Handling  - Failure  cases  are  flagged  by  setting  TAUAVG  = 10000, 
so  as  to  print  asterisks. 

9.  Flow  Diagram  - none 


10. 


Symbol  List  - none 


FKBOLT 


I 


t 


J 


1.  Algorithm  - contains  properties  for  bolts. 

2.  Argument  List  - KBOLT 

3.  Common  - DBLBLT 

4.  Length  - 50  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - DBLB,  STPLP 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List 

FT  - ultimate  tensile  strength 

FSS  - ultimate  shear  stress 

BK  - density  factor 

BG  - shear  modulus 


I 
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FKPROP 


1.  Algorithm  - Determine  the  material  properties  as  a function  of  W and  D. 

2.  Argument  List  - WD,  D 

3.  Common  - DBLBLT 

4.  Length  - 125  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - DBLB,  ANAL,  TH,  STPLP 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 
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FPROP 


1.  Algorithm  - Sets  up  graphite  material  properties  that  remain  fixed. 

2.  Argument  List  - MATL 

3.  Conmon  - DBLBLT 

4.  Length  - 52  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - DBLB,  ANAL,  STPLP 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List: 

PROP  (1)  = Ultimate  Tensile  Strength 

PROP  (2)  = Ultimate  Bearing  Strength 

PROP  (5)  = Elastic  Modulus 

MATL  * Material  Code 
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INIT 


1.  Algorithm  - Inltlatlzer  WORK  array  to  zero. 

2.  Argument  List  - none 

3.  Common  - WORK 

4.  Length  - 6 Octal 

5.  Subroutine  Called  - none 

7.  Subroutined  Called  By  - EDSEL,  PRINT,  B0LT4,  BOND5 

6.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 
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MAIN 


1.  Algorithm  - This  is  the  main  routine  that  initializes  data,  sets  the 
tablet  flag,  displays  the  three  main  options  and  branches  to  the  selected 
routine.  Upon  EXIT,  user  messages  are  displayed. 


i i 2.  Arguments  - none 

3.  Common  - blank,  CONSOL,  TABFLG,  BLK1 , BLK2 , WK 

4.  Length  - 417  Octal 

: I 

5.  Subroutines  Called  - OVERLAY  (DESIGN,  PRINT) 

6.  Subroutines  Called  By  - none 

[ i 

t j 7.  1/0  = If  data  exists  on  SAVE  file,  the  file  is  read  to  initialize 

\ the  ANAME  and  I TYPE  arrays. 

j Files  1 and  2 are  rewound. 

\ 

| The  main  option  list  is  displayed  on  unit  6 and  the  option  read. 

18.  Error  Handling  - An  Invalid  code  is  re-entered. 

If  the  user  states  that  data  exists  on  an  empty  file,  reading  that 
^ file  will  cause  a system  read  error  and  abort. 

| 

• 9.  Flow  Diagram  - see  Figure  15. 

F 

f 10.  Symbol  List  - The  following  are  not  covered  by  the  general  list: 

> ICODE  - main  option  code  entered 

\ ICON/CONSOL/  - set  = 1 If  consolidation  of  SAVE  file  selected 

| I TAB/TAB FLG/  - set  = 1 if  XYLOC  routine  to  use  tablet  calls 
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/ DATA  % 
ON  SAVE 
\ FILE 


WRITE 

H 

EXIT 

'•/RETURN  _''N' 

MESSAGES 

Figure  15. 

MAIN  Flow  Diagram 
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NAME 


1.  Algorithm  - This  routine  prompts  the  user  for,  and  reads  the  name 
assigned  to  each  analysis  problem,  and  displays  elapsed  time. 


i 


2.  Argument  List  - none 

3.  Common  - BLK1 

4.  Length  - 50  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - B0LT1 , B0LT4,  BONDI,  B0ND2,  B0ND5 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


L 


OPTION 


1.  Algorithm  - Draws  the  boxes  for  RETURN,  EXECUTE,  and  RE-INPUT, 
requests  a screen  coorlnate,  and  returns  a code  of  1 , 2 or  3,  respectively. 

2.  Argument  List  - IX,  IY,  IDEL 

3.  Common  - none 

4.  Length  - 174  Octal 


5.  Subroutines  Called  - XYLOC,  BOX 


6.  Subroutines  Called  By  - B0LT1 , B0LT4 


7.  1/0  - displays  boxes  with  labels  on  screen. 


8.  Error  Handling  - Keeps  requesting  a coordinate  until  one  is  within 
a box. 


9.  Flow  Diagram  - none 


Symbol 


IDELX 


delta  Y coordinate  of  box 

delta  X coordinate  from  left  side  of  box  for  labels 


IDEL2  - 2 * IDEL;  delta  X between  box  top  left  corners 

IR  - return  code 

IX  - starting  X coordinate  for  first  box 

I XX  - requested  X coordinate 

1X1  - starting  X coordinate  for  first  label 

1X2  - starting  X coordinate  for  second  label 

1X3  - starting  X coordinate  for  third  label 

IYY  - starting  Y coord,  for  first  box 

IY  - requested  Y coordinate 
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0UT1 


1.  Algorithm  - This  Is  the  output  routine  for  the  bolted  double-lap 
and  supported  single-lap  joints.  The  Input  argument  Is  the  unit  number 
for  the  formatted  writes. 


The  type  of  analysis  is  dependent  on  the  original  parameters  MO  and 
NX,  as  shown. 

2.  Argument  List  - NU 

3.  Common  - BLK1 , WORK 

4.  Length  - 452  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  BOLT!,  SELECT 

7.  1/0  - Writes  the  formatted  data  onto  either  the  display  (Unit  6), 
or  the  PRINT  file  (Unit  2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - See  Figure  16. 

10.  Symbol  List  - The  following  are  not  covered  by  either  the  general 
list  or  the  WORK  array  for  these  types: 

D6  - 6 * D 

FAIL  - array  of  possible  failure  modes 
IPB  - % equivalent  of  PB  decimal  array 

TYPE  - bolt  type  array 

WD  - W/D 
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! Figure  16.  0UT1  Flow  Diagram 

i 
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0UT4 

1.  Algorithm  - Output  routine  for  the  bolted  stepped-lap  joint. 

If  NU  = -6,  the  routine  has  been  called  to  view  only  existing  input 
data  contained  in  the  WORK  array. 

2.  Argument  List  - NU 

3.  Common  - BLK1 , WORK 

4.  Length  - 370  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  B0LT4,  SELECT 

7.  1/0  - Writes  the  formatted  data  to  either  the  display  (Unit  6),  or 
the  PRINT  file  (Unit  2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - The  following  are  not  covered  by  either  the  general 
list  or  the  WORK  array  for  bolted  stepped-lap  joints. 

BOLT  - bolt  type  array 

FAIL  - array  of  failure  modes 

INPUT  - flag  that  only  input  is  to  be  written 

NU  - unit  number 


5 

5 


j 
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0UT7 

1.  Algorithm  - Output  routine  for  solutions  of  bonded  double-lap  and 
supported  single-lap  joints. 

The  variable  NC  selects  the  appropriate  heading,  and  K the  type  of 
load  and  pertinent  properties. 

2.  Argument  List  - none 

3.  Common  - blank,  BLK1 , WORK 

4.  Length  - 676  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  BONDI,  SELECT 

7.  1/0  - Writes  the  formatted  data  to  either  the  display  (NU  = 6)  or 
the  PRINT  file  (NU  = 2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - See  the  appropriate  WORK  equivalence  descriptions. 
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0UT8 
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1.  Algorithm  - Output  routine  for  a bonded  unsupported  single-lap  joint. 
After  writing  out  the  Input  data,  the  output  formats  are  dependent  on 
the  Input  overlap  length,  W0RK(2).  If  zero,  seven  i/t  ratios  are  used, 
and  each  of  the  data  arrays  within  WORK  are  written  according  to  their 
respective  rormats. 

If  W0RK(1)  f 0,  strengths  are  written;  If  W0RK(1)  = 0,  stresses  are 
written. 

2.  Argument  List  - NU 

3.  Common  - BLK1 , WORK 

4.  Length  - 525  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  B0ND2,  SELECT 

7.  1/0  - Writes  the  formatted  data  to  either  the  display  (NU  = 6),  or 
PRINT  file  (NU  = 2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - Reference  the  general  list  for  variable  descriptions. 


74 


0UT9 

1.  Algorithm  - Output  routine  for  bonded  scarf  joints,  symnetrical 
and  asymmetrical. 

W0RK(1)  is  the  load  type  (-1,0,+!),  and  determines  the  value  of  K. 

W0RK(2)  determines  whether  strengths  or  stresses  were  computed. 

W0RK(3)  determines  whether  the  specified  overlap  or  a range  of  7 
overlaps  were  used. 

2.  Argument  List  - NU 

3.  Common  - BLK1 , WORK 

4.  Length  - 440  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  BONDI,  SELECT 

7.  1/0  - Writes  the  formatted  data  to  either  the  display  (NU  = 6) , or 
the  PRINT  file  (NU  = 2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List: 

END  - critical  end  array  (left,  right) 

K - load  type,  IWORK(l)  + 2 

NUM  - number  of  Items  In  each  array 

TYPE  - array  of  descriptive  load  types 
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0UT11 


1.  Algorithm  - Output  routine  for  bonded  stepped-! ap  joints  and  doublers. 
If  NU  is  negative,  only  the  input  data  is  to  be  displayed  for  the  user. 

The  compressed  format  and  variable  length  of  the  arrays  within  WORK 
is  due  to  the  variable  number  of  steps  for  the  different  computations. 
Reference  Table  2.  for  the  typical  layout  of  the  WORK  array. 

I W0RK( 15)  contains  number  of  segments  within  WORK,  from  1 to  4.  The 
first  is  the  input  step  data,  and  the  others  are  up  to  3 segments  of  com- 
putational output. 

2.  Argument  List  - NU 

3.  Common  - BLK1 , WORK 

4.  Length  - 626  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PRINT,  B0ND5,  SELECT 

7.  I/p  = Writes  formatted  data  to  either  the  display  (NU  = 6),  or  the 
PRINT  file  (NU  = 2). 

8.  Error  Handling  - none 

9.  Flow  Diagram  - Figure  17 
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10.  Symbol  List  - The  following  describes  the  unique  variables  used. 

Refer  to  the  general  list  and  the  WORK  equivalence  descriptions  for 

i 

the  remal nder : 

CST  - array  of  load  types 

IL  - Initial  location  In  each  WORK  segment 

IT  - load  type  pointer  for  CST 

IUPDAT  - set  = 1 If  only  Input  to  be  displayed 

KA  - no.  of  arrays  within  WORK  for  previous  segment 

M - no.  of  segments 

N - no.  of  steps  analyzed  for  each  segment 

i 

t 


PCT 


1.  Algorithm  - This  routine  determines  the  bolt  load  distribution  for 
up  to  6 rows  of  bolts  of  a double-lap  joint. 


2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 


Argument  List  - M,  WD,  D,  T,  PB,  A 
Common  - DBLBLT 
Length  - 317  Octal 
Subroutines  Called  - none 

Subroutines  Called  By  - DBLB,  ANAL,  WDMAX,  WDMIN,  TH 
1/0  - none 

Error  Handling  - none 
Flow  Diagram  - none 
Symbol  List  - none 


'i 

i 

jj 
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PCSTP 


[ 

i 

! 
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1.  Algorithm  - This  routine  calculates  the  amount  of  load  transferred 
by  each  row  of  bolts  In  a stepped-lap  joint. 

2.  Argument  List  - El,  EO,  G,  M,  S,  D,  W,  TI,  T0,  PCT 

3.  Common  - none 

4.  Length  - 452  Octal 

5.  Subroutines  Called  - SID 

6.  Subroutines  Called  By  - STPLP 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


PRINT 


'I 


[ i 1.  Algorithm  - This  routine,  so  named  because  it  originally  contained 

only  the  option  to  print,  also  consolidates  the  SAVE  file. 

If  the  user  has  selected  the  selective  output  option  (ICON  =0),  he 
is  asked  whether  the  solution  is  to  be  displayed  or  output  to  the  PRINT 
file. 

After  the  appropriate  title,  all  the  names  on  the  SAVE  file,  and 
option  boxes  displayed,  the  user  is  prompted  to  pick  the  names  desired 
by  their  screen  location. 

Execution  of  the  print  option  is  accomplished  by  testing  the  IPRNT 
array,  locating  the  solution  on  the  SAVE  file,  reading  the  data  into  WORK, 
and  calling  the  appropriate  output  routine.  If  the  solution  is  displayed, 
XYLOC  is  called  to  interrupt  the  loop  until  the  user  is  ready  to  continue. 

Execution  of  the  consolidation  option  consists  of  copying  all 
solutions,  except  those  flagged,  from  the  SAVE  file  to  a temporary  file 
(TAPE3)  then  back  to  the  SAVE  file.  The  design  numbers  are  resequenced 
in  the  process. 

2.  Argument  List  - none 

3.  Comrton  - blank,  BLK1 , BLK2,  CONSOL,  WORK 

4.  Length  - 1251  Octal 

5.  Subroutines  Called  - INIT,  0UT1,  0UT4,  0UT7,  0UT8,  0UT9,  0UT11 , 

READ1 , XYLOC 
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6.  Subroutine  Called  By  - MAIN 


1 

7.  1/0  - Units  5 and  6 are  used  to  read  and  display  data.  For  the  con- 
solidation option,  the  selected  contents  of  Unit  1 are  copies  to  Unit  3; 
Unit  3 is  then  copied  back  to  Unit  1,  to  complete  consolidation  of  the 
SAVE  file. 

8.  Error  Handling  - If  EXECUTE  has  been  selected  without  any  names 
picked,  the  user  is  requested  to  make  another  coordinate  selection. 


9.  Flow  Diagram  - Figure  18 

10.  Symbol  List: 

IPRNT 

- array  for  flagged  solutions 

IT 

- solution  type 

IXMIN 

- min.  X coord,  of  a name  on  screen 

IXP 

- array  of  starting  X locations  for  screen  names 

IYMAX 

- max.  (top)  Y coordinate  for  name  detects 

IYMIN 

- min.  (bottom)  Y coordinate  for  name  detects 

IYP 

- array  of  Y coordinate  for  up  to  10  rows  of  names 

KC 

- column  number  of  selected  name  (from  1-10) 

KD 

- row  number  of  selected  name  (from  1-10) 

KN 

- design  number  selected  (10*KD+KC) 

KNEW 

- counter  for  NOT  during  consolidation 

: 

KOUT 

- output  response  for  display 

LOOP 

- flag  for  consolidation 

NA 

- colsolidation  read  unit 

NB 

- consolidation  write  unit 

NU 

- print  unit  number 
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Figure  18.  PRINT  Flow  Diagram 
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Figure  18  PRINT  Flow  Diagram  (continued) 
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PSCARF 


1 

1.  Algorithm  - This  routine  Is  for  the  elastic-plastic  analysis  of 
unbalanced  scarf  joints. 

2.  Argument  List  - OL,  ETR,  THERMC , GAMMAR,  LCRTND,  TAUAVG 

3.  Common  - none 

4.  Length  - 1105  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - ANAL9 

j 7.  1/0  - none 

i 

8.  Error  Handling  - If  TAUAVG  > 1,  there  Is  an  error;  TAUAVG  is  set 
equal  to  1000  as  a flag  and  control  returned  to  the  calling  routine. 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 

f 


r 
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QUADMN 


f. 


1.  Algorithm  - This  routine  finds  the  X associated  with  the  minimum 
f(x)  by  quadratic  Interpolation. 

2.  Argument  List  - X,  F,  OX,  XMAX,  DDX,  IERR 

3.  Common  - none 

4.  Length  - 525  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - DBLB 

7.  1/0  - none 

8.  Error  Handling  - The  following  error  codes  may  be  returned  by  IERR. 

0 = no  error 

1 = XO  > (XMAX-DX)  and  DFO  > 0. 

2 = XO  > XMAX 

3 = iterations  exceeded  200  max. 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 
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READ1 

1.  Algorithm  - This  routine  searches  the  SAVE  file  for  the  requested 
design  number.  If  greater  than  100,  the  file  is  read  to  the  end  of  data. 
If  less  than  100,  the  routine  locates  the  solution,  reads  the  name  and 
type,  and  returns. 

2.  Argument  List  - IDES 

3.  Common  - blank,  BLK1 , BLK2 

4.  Length  - 70  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - EDSEL,  SAVE,  PRINT 

7.  1/0  - Reads  Unit  1 for  desired  design  number. 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List: 

IDES  - requested  analysis  design  number 
KD  - design  number  read  from  SAVE  file 
KT  - design  type  read  from  SAVE  file 
TEMP  - design  name  read  from  SAVE  file. 


SAVE 


1.  Algorithm  - When  a solution  Is  selected  for  saving  on  the  SAVE 
file  (Unit  1),  this  routine  Is  called  . It  only  allows  100  solutions  to 
be  output.  After  writing  the  WORK  array  data  to  the  SAVE  file,  the 
end-of-data  record  is  written  and  a backspace  done  to  force  a buffer  dump, 
and  position  the  file  for  the  next  save. 

2.  Argument  List  - IT,  NENT 

3.  Comnon  - blank,  BLK1 , BLK2,  WORK 

4.  Length  - 142  Octal 

5.  Subroutines  Called  - READ! 

6.  Subroutines  Called  By  - SELECT  I 

j 

7.  1/0  - Writes  2 records  to  the  SAVE  file  containing  the  design  number,  j 

analysis  name,  analysis  type,  and  the  work  array  entries.  A third  record  \ 

j 

is  then  written  containing  the  end-of-data  flag.  j 

j 

j 

8.  Error  Handling  - If  the  number  of  solutions  equals  100,  a message  is  j 

displayed  that  the  save  was  aborted  due  to  max.  solutions. 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 

i 
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SELECT 

1.  Algorithm  - This  routine  processes  the  user  options  after  display 
of  an  executed  analysis. 

The  options  to  PRINT,  SAVE,  re-analyze,  and  return  are  displayed  on 
a line,  separated  by  asterisks.  The  user  Is  requested  to  make  selections 
by  the  horizontal  screen  location  picked.  The  PRINT  and  SAVE  options  must 
be  processed  before  a re-analysis  or  return  option  for  obvious  reasons. 
When  either  of  the  first  two  options  Is  processed,  a flag  Is  set  to 
prevent  repeating  that  selection  Inadvertently. 

If  RETURN  Is  selected,  the  value  of  N Is  returned  as  zero  to  the 
calling  routine. 

2.  Argument  List  - N,  NENT 

3.  Common  - un labeled 

4.  Length  - 174  Octal 

5.  Subroutines  Called  - 0UT1 , 0UT11,  0UT4,  0UT7,  0UT8,  0UT9,  SAVE,  XYLOC 

6.  Subroutines  Called  By  - B0LT1 , B0LT4,  BONDI,  B0ND2,  B0ND5 

7.  1/0  - display  of  options 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List 

IP  - set  = 1 after  PRINT  selection 

IS  - set  = 1 after  SAVE  selection 

N - Input  = analysis  type 

Output  = return  code  (0  = RETURN) 
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1.  Algorithm  - A single-precision  simultaneous  equation  solver, 
inverse  finder,  and  determinent  routine. 

2.  Argument  List  - A,  N,  NDROW , NDCOLA , B,  M,  NDCOLB,  SIGDIG,  I ERROR, 

PIVOT,  INDEX,  SCALEB 

3.  Common  - none 

4.  Length  - 607  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - PCTSTP 

7.  1/0  - none 

8.  Error  Handling  - IERR0R  returned  to  calling  routine 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 
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STPLP 

1.  Algorithm  - Bolted  stepped-lap  joint  analysis  routine  to  determine 
the  amount  of  load  retained  by  the  Inner  and  outer  adherends,  and  determine 
the  margin  of  safety  and  failure  mode. 

i 2.  Argument  List  - none 

3.  Common  - WORK,  BLOCK 

4.  Length  - 377  Octal 

5.  Subroutines  Called  - FKPROP,  FKBOLT,  FPROP,  PCTSTP 

6.  Subroutines  Called  By  - B0LT4 

i 

i 7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 
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— •"  - 1 

TO 

Algorithm  - Function  subprogram  for  determining  the  joint  thickness 

\ 

• • of  double-lap  bolted  joint. 

f 2. 

Argument  List  - NX,  M,  WD,  D 

\ 

I 3. 

Common  - DBLBLT 

j 4. 

Length  - 340  Octal 

, 5‘ 

Subroutines  Called  - FKPROP,  PCT 

f 

6. 

Subroutines  Called  By  - DBLB,  WDMAX,  WDMIN,  WWT 

7. 

1/0  - none 

i 

8. 

Error  Handling  - none 

i 9. 

Flow  Diagram  - none 

► . 

10. 

Symbol  List  - none 

■ 

k. 1 

: 

92  1 

k 

WDMAX 


r 


f 

i 

! 


I 

E 

I 

[ 


i 

f 

t 

t 

l 

l 


1.  Algorithm  - Determines  a maximum  W/D  ratio  for  a bolted  double-lap 
joint. 

2.  Argument  List  - NX,  M,  WD,  D 

3.  Coirmon  - DBLBLT 

4.  Length  - 243  Octal 

5.  Subroutines  Called  - PCT,  TH 

6.  Subroutines  Called  By  - DBLB 

7.  1/0  - None 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


\ 
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WDMIN 

1.  Algorithm  - Determines  a minimum  W/D  ratio  for  a bolted  double-lap 
joint. 

2.  Argument  List  - NX,  M,  WD,  D 

3.  Coirmon  - DBLBLT 

4.  Length  - 265  Octal 

5.  Subroutines  Called  - PCT,  TH 

6.  Subroutines  Called  By  - DBLB 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


WT 


1.  Algorithm  - Function  subprogram  that  determines  the  weight  penalty 
for  a double-lap  bolted  splice. 


i 

! 


2.  Argument  List  - X 

3.  Common  - OBLBLT 

4.  Length  - 66  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - DBLB 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


I 
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WTSL 


Algorithm  - Determines  the  weight  penalties  for  single-lap 
bolted  splice. 

Argument  List  - X 

Common  - DBLBLT 

Length  - 44  Octal 

Subroutines  Called  - none 

Subroutines  Called  By  - DBLB,  WWT 

1/0  - none 

Error  Handling  - none 
Flow  Diagram  - none 
Symbol  List  - none 
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WWT 


1.  Algorithm  - A function  that  sets  up  the  X array  for  determining  the 
weight  penalty  of  a bolted  splice. 

2.  Argument  List  - WD 

3.  Common  - OBLBLT 

4.  Length  - 43  Octal 

5.  Subroutines  Called  - WT,  TH  , WTSL 

6.  Subroutines  Called  By  - DBLB 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - none 


XYLOC 

1.  Algorithm  - This  routine  calls  PL0T10  routines  for  screen  X & Y 
coordinates.  If  ITAB  = 0,  the  routine  calls  DCURSR  for  the  coordinates. 

If  not,  the  user  is  notified  that  the  tablet  Is  on;  the  routine  then  loops 
to  track  the  cursor  position  on  the  screen  until  the  pen  Is  depressed 

on  the  tablet. 

2.  Argument  List  - IX,  IY 

3.  Common  - TABFLG 

4.  Length  - 76  Octal 

5.  Subroutines  Called  - none 

6.  Subroutines  Called  By  - SELECT,  EDSEL,  OPTION,  PRINT 

7.  1/0  - none 

8.  Error  Handling  - none 

9.  Flow  Diagram  - none 

10.  Symbol  List  - 

IC  - dummy  ASCII  character 
IH  - ASCII  character 
IX  - X screen  coord. 

IY  - Y screen  coord. 


BOATED  — STANDARD  DOUBLE  LAP  JOINT 
ANALYSIS  NAME  - B0LT1-9 
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Figure  19.  Bolted  Oouble-Lap  Input  (Analysis) 


Figure  20.  Bolted  Double-Lap  Output  (Analysis) 
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Figure  22.  Bolted  Double-Lap  Output  (Optimization) 
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Figure  24.  Bolted  Double-Lap  Output 


BOLTED  — STANDARD  DOUBLE  LAP  JOINT 
ANALYSIS  NAME  = BOLT1-13 


Figure  25.  Bolted  Double-Lap  Input  (P/0,  N?*0) 
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Figure  26.  Bolted  Double-Lap  Output  (Margins  and  Weight) 
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Figure  27.  Bolted  Unsupported  Single-Lap  Input 
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Figure  28.  Bolted  Unsupported  Single-Lap  Output 


BOLTED  — SUPPORTED  SINGLE-LAP  JOINT 
ANALYSIS  NAME  = BOLT-3 


Figure  29.  Bolted  Supported  Single-Lap  Input 
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Figure  30.  Bolted  Supported  Single-Lap  Output 
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Figure  31.  Bolted  Stepped-Lap  Input 
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Figure  32.  Bolted  Stepped-Lap  Editing 


fOt TCP  — STtfftS-L*»  JO  JUT  ANALYSIS 
ANALYSIS  NANE  • S0LT4-J4 


®HIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  COPY  FURNISHED  TO  DDC  _ 


Figure  33.  Bolted  Stepped-Lap  Output 
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Figure  34.  Bolted  Stepped-Lap  Modifications 


•OLTEO  — ST£NP£D-U*N  JOINT  W1ALTSIS 
OdOLVSIS  NONE  • E01T4-D 


1 

I 

! 


, 


> 


i- 


i i 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURNISHED  TO  DDC 


IS 

tt 

IS 


HM 

tt 

jS3 

AM 


• 88 


-v, 

Xkh  UJUI 
« O MM 
JN  •«  JJ 

“ *Ok  Sl a. 


21^5  jyy 


J*j6z 

---Sstttt 

->oauta33 

2855*8* 


y--*fU 

OKhh 

r r z k mm 


Vise 

Wki 

Is »! 


H if 


J V c-ni 
8*  ** 


J 

2* 


if  £2 


$ Is 


: &5 


£ 


tt 

fun 


gt! 

vi22 


•j** 

>oc* 


st 

aim 


0>> 

&■  • 

Coo 


2 

* 


S 


15 


ecu 


JJJC 

8882 

MNEE 
I • • I 


| fi* 

M 

2 2f 


J?E6  . 

2 = Sc 


8 <•! 


115 


tt 

*• 


w SS 

m mm 

£ SS 
tttt 


1 

tt 


s SI  - 


s 

t 

g 


J 


Figure  35.  Bolted  Stepped-Lap  Output 
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Figure  36.  Bonded  Double-Lap  Selection 


BONDED  — STANDARD  DOUBLE-LAP  JOINT 
ANALYSIS  NAME  * BONDS-15 


Figure  37.  Bonded  Double-Lap  Input  (All) 
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Figure  38.  Bonded  Double-Lap  Output 


Figure  40.  Bonded  Double-Lap  Re-Display  & Modify 
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Figure  41.  Bonded  Double-Lap  Output 
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Figure  42.  Bonded  Unsupported  Single-Lap  Input 


IBIS  PAGE  IS  BEST  QUALITY  PEACTICABlil 
EBQM  OGPY  ElWwliiHJiJJ  TO  DUG  - 

5 


TED  SINGLE-LAP  JOINT  ANALYSIS 
ALVSIS  NANE  • ION01S-0 


Figure  45.  Bonded  Unsupported  Single-Lap  Output  (P=0,  01=0) 
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Figure  46.  Bonded  Unsupported  Single-Lap  Output  (P/0,  0L=0) 
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Figure  47.  Bonded  Supported  Single-Lap  Input 


Figure  48.  Bonded  Supported  Single-Lap  Output 
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Figure  49.  Bonded  Stepped-Lap  Joint  Input 
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Figure  50.  Bonded  Stepped-Lap  Joint  Re-Input  (Editing) 
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Figure  51.  Bonded  Stepped-Lap  Joint  Output  - Input  Data 
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Figure  52.  Bonded  Stepped-Lap  Joint  - Elastic  Analysis 
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Figure  53.  Bonded  Stepped-Lap  Joint  - Elastic-Plastic  Analysis 
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Figure  54.  Bonded  Stepped-Lap  Joint  Output  - El. -PI.  Analysis  With  Infinite  Adherend  Strength 
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Figure  55.  Bonded  Stepped-Lap  Joint  Output  Summary 
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Figure  56.  Bonded  Scarf  Joint  Input 
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Figure  57.  Bonded  Scarf  Joint  Output  (Load^O,  Overlap^) 
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Figure  58.  Bonded  Scarf  Joint  Output  (Load=0,  0verlap=0) 
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Figure  59.  Bonded  Scarf  Joint  Output  (Load*0,  OverlapfO) 
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Figure  60.  Selective  Output  Processing  Option 


SELECTIVE  OUTPUT  OF  SOLUTIONS  TO  PRINT  FILE 


Figure  61.  Selective  Output  of  Solutions  To  Print  File 
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Figure  62.  Selection  of  Save  File  Solutions  For  Display 


CONSOLIDATION  OF  SAUE  FILE  SOLUTIONS 


Figure  64.  Consolidate  Solutions  On  Save  File 


SECTION  II 


ANALYTICAL  PROBLEM  DESCRIPTION  FOR  BONDED  JOINTS 

PHYSICAL  PROBLEM  DESCRIPTION 

Figure  65  Illustrates  the  various  types  of  bonded  joint  geometries 
covered  by  this  program.  The  difference  between  a joint  and  a doubler  is 
that,  in  a joint,  the  entire  load  is  passed  through  the  bond  whereas  in 
a doubler,  only  some  of  the  load  Is  so  transferred. 

The  loads  in  the  adherends  are  basically  in-plane,  with  no  lateral  applied 
loads.  However,  in  the  overlap  area,  there  can  be  bending  deflections  and 
transverse  stresses  in  one  or  more  of  the  adherends  as  a result  of  eccentricities 
in  load  path.  The  analyses  cover  remote  loads  which  are  tensile,  compressive, 
or  In-plane  shear. 

The  adherends  are  treated  as  homogeneous  orthrotropic  materials  rather 
than  as  multi -component  fiber/matrix  combinations.  Experimentally  determined 
material  stress-strain  curves  for  the  adhesive  and  adherends  are  used  to  cover 
the  possible  failure  modes  of  adhesive  shear,  adhesive  peel,  adherend  in-plane 
failure  away  from  the  bonded  overlap,  and  adherend  interlaminar  failure  in 
the  joint  area.  Thermal  properties,  as  well  as  mechanical,  are  provided  for 
so  that  the  analyses  can  account  for  the  initial  stresses  induced  by  the 
bonding  of  metals  to  composites.  The  adherends  are  treated  as  linearly 
elastic  to  failure,  while  the  adhesive  Is  modelled  as  an  elastic-plastic 
material  to  account  for  its  considerable  non-linear  deformation  prior  to 
failure.  The  linear  treatment  of  the  adherends  does  not  Impose  any  real 
restriction  on  the  utility  of  the  programs  because  a sustained  yield  load  on 
the  adherends  causes  steady  and  progressive  failure  of  the  bond  as  long  as 
that  load  is  maintained. 
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DOUBLE -STRAP  JOINT 


MATHEMATICAL  MODEL  DESCRIPTION 


t \ 1 O 

f The  mathematical  models  used  are  fully  described  , and  the  features 

accounted  for  here  are  outlined  in  the  physical  description.  No  published 
solutions  cover  all  the  governing  variables  precisely  but  it  is  felt  that 
those  of  these  analyses  are  the  most  important.  The  importance  of  accounting 
for  adhesive  non-linearity  is  paramount. 

One  key  improvement  found  in  these  solutions  and  not  found  elsewhere  is 
the  accounting  for  adhesive  plasticity  in  an  explicit  simple  manner.  Other 
more  elaborate  characterizations  have  been  tried  elsewhere.  This  has  pre 
; vented  some  solutions  from  being  completed  to  workable  form  and  made  others 

too  timecconsuming  for  production  work.  An  explicit  closed  form  analysis  of 
double-lap  bonded  joints1  has  shown  that  any  bi-linear  (two  straight  lines) 
representation  of  adhesive  shear  characteristics  leads  to  precisely  the  same 
j joint  strength  provides  that  the  adhesive  models  have  the  same  failure  stress, 

! strain,  and  strain  energy.  Only  the  elastic-plastic  limit  of  this  family  of 

s i 

solutions  has  led  to  explicit  closed-form  solutions  for  more  complex  joint 
geometries  and  this  is  why  it  has  been  adopted  for  this  work. 

; A second  key  improvement  pertains  only  to  the  unsupported  single-lap 

joints.  This  is  in  the  calculation  of  the  bending  moment  at  the  ends  of  the 

overlap.  Since  this  serves  as  a dominant  boundary  condition  on  the  adhesive 

shear  stresses  and  peel  stresses  as  well  as  the  adherend  bending  stresses, 

it  must  be  determined  as  accurately  as  possible.  With  the  exception  of  one 

2 

Ph.D.  thesis  which  handles  this  detail  to  the  same  accuracy  , all  other 
published  solutions  for  single-lap  bonded  joints  ase  unnecessarily  coarse 
| approximations  in  determining  that  bending  moment.  This  nullifies  the 


1.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double-Lap  Joints",  NASA  CR-112235, 

■ January  1973. 

2.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Single-Lap  Joints",  NASA  CR-1 12236, 

' January  1973. 

3.  Hart- Smith,  L.  J.,  "Adhesive-Bonded  Scarf  and  Stepped-Lap  Joints", 

NASA  CR-1 12237,  January  1977. 
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fc 

r 
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quantitative  value  of  the  subsequent  portions  of  those  investigations,  even 
when  they  are  qualitatively  correct  elsewhere. 

Other  published  bonded  joint  solutions  contain  up  to  two  factors  not 
accounted  for  here.  One  is  the  ply-by-ply  stress  distribution  within  the 
laminate,  as  opposed  to  the  homogeneous  orthotropic  model  used  here.  Published 
assessments  of  this  effect  indicate  that  provision  for  interlaminar  shear 
deformation  alone  does  not  have  a major  impact  on  the  adhesive  shear  stress 
distribution  but  that  the  variation  of  the  interlaminar  tension  stresses 
through  the  thickness  does  have  a significant  effect  on  the  elastic  adhesive 
shear  stresses.  The  result  of  accounting  for  this  factor  is  that  the  adhesive 
shear  stresses  peak  slightly  inside  the  end  of  the  overlap,  instead  of  right 
at  the  end  where  equilibrium  demands  that  there  be  no  stresses.  This  effect 
is  particularly  significant  for  unsupported  single-lap  joints  and  is  not 
negligible  for  double-lap  joints  either,  provided  that  attention  is  restricted 
to  linearly  elastic  adhesive  behavior.  Near  the  ultimate  joint  strength, 
however,  the  adhesive  non-linear  behavior  makes  this  effect  less  and  less 
significant.  Certainly,  it  is  more  important  to  account  for  adhesive  plasticity, 
then,  than  for  the  variation  in  peel  stresses  through  the  thickness.  Account- 
ing for  the  latter  raises  the  order  and  complexity  of  the  governing  differen- 
tial equations.  So  far,  no  analyses  have  been  able  to  account  for  both  of 
these  effects  simultaneously. 

The  second  factor  treated  more  precisely  elsewhere  is  that  the  adhesive 
failure  criteria  here  are  for  shear  and  peel  separately,  instead  of  in  terms 
of  an  interaction  formula.  The  justification  for  the  present  simpler  approach 
is  as  follows.  When  peel  stresses  are  significant  in  comparison  with  the 
shear  stresses,  the  former  enforce  a joint  strength  reduction.  Therefore, 
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the  Interest  In  peel  stresses  Is  mainly  to  Identify  those  geometries  which 
need  modifying  to  alleviate  such  stresses.  In  the  satisfactory  joint  geometries, 
the  peel  stresses  can  be  neglected  in  comparison  with  the  shear  stresses  and, 
in  identifying  those  unsatisfactory  geometries  in  which  the  peel  stresses  are 
predominant,  the  shear  stresses  can  be  ignored.  There  is  only  a small  range 
of  geometries  in  which  both  stresses  are  significant  simultaneously  and  even 
those  can  be  improved  to  reduce  the  strength  loss  due  to  the  peel  stresses. 

The  solutions  are  formulated  in  terms  of  differential  equations  from 
classical  continium  mechanics.  Some  explicit  solutions  are  derived.  OthPrs 
are  exact,  but  implicit,  and  require  an  iterative  solution.  The  scarf  joint 
solutions  use  a finite  number  of  coefficients  of  a power  series  solution. 
Consequently,  the  scarf  joint  strengths  are  obtained  as  the  integral  of  these 
stresses,  with  considerable  accuracy,  but  the  internal  stress  distributions 
are  usually  not  sufficiently  accurate  and  are,  therefore,  not  specified.  The 
stepped-lap  programs  encounter  potential  numerical  accuracy  problems  because 
of  the  very  high  shear  stress  gradients,  which  sometimes  caused  failure  to 
converge.  All  such  cases  known  have  been  eliminated  by  artifically  dividing 
the  step  lengths  automatically  to  restrict  the  arguments  of  the  exponential 
functions.  The  internal  operations  are  tested  to  16  significant  digits  and 

the  overall  solutions  are  accurate  to  at  least  6 significant  digits. 

DESCRIPTION  OF  NUMERICAL  METHODS 

1-3 

The  details  of  the  analysis  methods  are  fully  documented. 


1.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double-Lap  Joints",  NASA  CR-112235, 
January  1973. 

2.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double-Lap  Joints",  NASA  CR-11236, 
January  1973. 

3.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Scarf  and  Stepped-Lap  Joints", 

NASA  CR- 112237,  January  1977. 
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LIMITATIONS 

Other  than  the  FORMAT  limitations,  there  are  no  mathematical  limitations 
on  these  programs.  The  input/output  Instructions  are  set  up  for  U.S.  customary 
units,  rather  than  S.I.  units,  and  changes  would  be  necessary  to  accommodate 
other  units. 

SOLUTION  ACCURACY 

The  solution  printed  will  be  accurate  to  at  least  the  number  of  signifi- 
cant digits  printed  with  only  one  exception.  That  exception  is  the  case  where 
the  adherends  are  so  thick  and  the  thermal  mismatch  so  great  that  the  joint 
breaks  apart  without  the  application  of  any  mechanical  load.  In  most  such 
cases,  this  situiation  is  identified  by  obviously  self-inconsistent  answers 
and  by  abnormally  long  run  times  because  of  failure  to  converge.  A concealed 
failure  case  is  that  in  which  the  stiffness  imbalance  and  thermal  mismatch 
cancel  each  other  out  for  one  particular  load  direction.  In  some  such 
instances,  the  removal  of  the  mechanical  load  result  in  failure  due  to  the 
thermal ly-induced  stresses.  There  is  no  special  provision  involved  automati- 
cally for  such  a case.  However,  re-running  the  problem  for  a near  zero  load 
or  reversed  temperature  differential  (to  simulate  reversal  of  the  load)  will 
bring  such  situations  to  light. 

DEFINITION  OF  NOTATION 

1-3 

The  notation  used  is  fully  explained. 


1.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double-Lap  Joints",  NASA  CR-112235, 

January  1973. 

2.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double-Lap  Joints",  NASA  CR-11236, 

January  1973.  ; 

3.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Scarf  and  Stepped-Lap  Joints", 

NASA  CR-112237 , January  1977. 
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ANALYTICAL  PROBLEM  DESCRIPTION  FOR  BOLTED  JOINTS 

PHYSICAL  PROBLEM  DESCRIPTION 

Four  types  of  bolted  joints  are  treated  here.  These  joints  represent 
some  of  the  more  common  techniques  of  joining  major  load  carrying  structure 
constructed  with  an  advanced  composite  material  of  graphite-epoxy.  These 
four  joint  types  which  are  illustrated  in  Figure  66  are: 

o balanced  double-lap 
o supported  single-lap 
o unsupported  single-lap 
o stepped- lap 

The  double-lap  joint  is  one  of  the  most  common  of  all  bolted  structural 
joints.  This  type  of  joint  is  well  suited  for  a wing  centerline  splice 
which  is  not  exposed  to  the  airstream.  It  is  a fairly  efficient  joint 
which  is  free  of  eccentrically  induced  bending  stresses.  It  is  also  a 
relatively  easy  and  inexpensive  joint  to  fabricate.  The  balanced  double- 
lap joint  is  essentially  a butted  joint,  symmetrically  spliced  by  two 
plates  of  equal  thickness.  Away  from  the  splice,  a nominal  skin  thickness 
Is  necessary  to  carry  the  design  tensile  load,  but  at  the  joint  this 
nominal  thickness  is  built  up  to  a thickness  which  can  carry  the  concen- 
trated stresses  at  the  bolts.  The  built-up  skins  are  then  butted  together 
and  spliced  by  two  Identical  plates,  each  of  half  the  thickness  of  the 
built-up  skin. 
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The  stepped-lap  joint  is  also  a joint  where  the  nominal  thickness  is 
built  up,  and  which  is  free  of  eccentrically  induced  bending  stresses.  But 
rather  than  having  a uniform  buildup  spliced  by  a uniform  splice  plate,  the 
splice  plate  is  stepped  and  fitted  to  a matching  stepped  buildup  as  shown 
in  Figure  66.  The  stepped-lap  joint  is  a more  efficient  joint  than  the 
double-lap  joint,  but  is  considerably  more  difficult  and  costly  to  fabricate 
and  fit.  It  is  important  to  appreciate  the  reason  for  the  greater  efficiency 
of  the  stepped-lap  joint.  First  of  all,  in  the  case  of  several  bolt  rows, 
with  each  row  on  a step,  the  transfer  load  distribution  is  more  uniform  than 
for  the  double-lap  joint  where  the  end  bolts  are  so  highly  loaded.  Secondly, 
and  more  significantly,  the  total  thickness  of  the  stepped-lap  joint  can  be 
less  than  the  double-lap  joint.  Where  all  of  the  load  is  in  the  buildup  and 
it  must  be  thickest,  the  splice  plate  need  only  be  minimal,  and  conversely, 
when  all  of  the  load  has  been  transferred  to  the  splice  plate,  the  buildup  need 
only  be  minimal.  As  a further  consequence,  the  bolt  lengths  can  be  shorter 
for  the  stepped-lap  joint. 

v 

Single-lap  joints  may  be  either  supported  or  unsupported.  No  bending 
stresses  are  assumed  to  exist  in  the  supported  single-lap  joint,  so  that 
the  structural  behavior  of  that  joint  is  nearly  identical  to  a double-lap 
joint  of  half  the  thickness  and  no  buildup.  The  support  for  this  type  joint 
is  provided  in  the  overlap  region  by  connecting  structure.  The  connecting 
structure  can,  for  example,  be  rib  flanges  supporting  the  overlap  region  of  a 
single-lap  skin  splice.  Eccentrically  induced  bending  stresses,  do,  however, 
exist  in  the  unsupported  single-lap  joint,  and  these  bending  stresses  are 
greatest  at  the  two  points  just  outside  the  overlap.  Single-lap  joints  are 
the  least  expensive  and  easiest  of  all  joint  types  to  fabricate.  Consequently, 
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the  added  expense  and  complexity  of  providing  a joint  buildup  is  generally 
rejected  in  favor  of  the  inefficient  uniform  thickness.  Of  course,  this 
thickness  is  considerably  more  than  the  nominal  skin  thickness  discussed 
previously.  Due  to  the  extreme  inefficiency  of  unsupported  single-lap  joints, 
their  applicability  should  be  reserved  for  low  tensile  load  carrying  secondary 

} 

structure  such  as  fairings,  or  perhaps  rib  and  spar  webs.  f 

MATHEMATICAL  PROBLEM  DESCRIPTION 

Development  of  a program  which  optimumly  designs  three  of  the  four  joint 
types  has  been  completed.  For  a given  joint  loading,  the  number  of  bolt 
rows,  bolt  diameter,  bolt  spacing,  and  joint  thickness  are  computed  for  the 
least  weight  joint.  This  search  for  the  optimum  design  may  be  conducted  for 
joints  of  the  following  type: 

o double-lap 

o supported  single-lap 

o unsupported  single-lap 

For  the  stepped-lap  joint  an  analysis  is  conducted  which  provides  margins  of 
safety  for  each  bolt  row  with  the  associated  failure  mode  for  a given  joint 
loading.  The  stepped-lap  joint  consists  of  too  many  variables  for  the  search 
algorithm  employed  for  the  other  three  joint  types. 


In  order  to  simplify  the  search  for  an  optimum  design,  the  number  of  independent 
variables  must  be  kept  to  a minimum.  For  this  reason,  the  math  model  associated 
with  the  three  types  of  joints  which  considered  here  for  optimum  joint  designs 
have  the  following  restrictions. 
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o same  bolt  sizes  for  all  rows 
o same  bolt  spacing  for  all  rows 
o edge  distance  of  3D 

o row  pitch  of  6D 

o two  compositions  of  graphite-epoxy  material 

o two  bolt  materials 

In  addition  the  double-lap  joint  is  restricted  to  a balanced  double-lap 
which  requires  that  the  two  outer  splice  plates  each  be  one-half  the  joint 
buildup  thickness.  Furthermore,  the  joint  buildup  is  assumed  to  be  gradual, 
with  a slope  of  15:1 . 

The  two  compositions  of  graphite-epoxy  material  were  chosen  to  differ  in 
the  percentage  of  0°  plies  since  the  stress  concentration  relief  factors 
for  these  composite  materials  were  found  to  be  predominantly  dependent  on  this 
percentage.  The  two  percentages  selected  were  25%  and  37-1/2%  of  the  0°  plies. 
The  25%  composition  of  course  corresponds  to  the  popular  pseudo-isotropic 
layup.  To  date,  stress  concentration  relief  data  has  been  accuimiulated  for 
these  two  compositions  only.  Of  course,  the  math  model  can  be  expanded  to 
include  more  material  compositions  whenever  this  data  becomes  available. 

The  stepped-lap  joint  analysis  allows  for  varying  step  lengths.  Bolt 
diameters  and  bolt  spacing  may  vary  from  one  step  to  another.  Step  thicknesses 
of  the  splice  plate  and  the  joint  buildup  are  independent  (i.e.,  the  joint 
need  not  be  balanced).  The  following  restrictions  exist  for  the  stepped-lap 
joint  math  model. 
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o bolts  at  center  of  steps 


o two  compositions  of  graphite-epoxy  material 

o two  bolt  materials 

Strictly  speaking,  the  bolts  of  all  the  rows  are  in  line  for  the  math 
models  employed  for  the  bolted  joints.  The  coupling  effects  of  the  bolt 
rows  has  been  neglected,  so  that  the  bolts  in  one  row  do  not  influence  the 
stress  state  at  another  row.  Therefore,  the  bolted  joint  analyses  developed 
here  are  not  strictly  valid  for  stepped-bolts,  but,  depending  on  the  row 
pitch,  may  not  incur  unreasonable  error. 

MATHEMATICAL  DERIVATION  OF  SOLUTIONS 

Four  failure  modes  are  treated  for  all  of  the  bolted  joint  designs  at 
each  bolt  row. 

o tension  at  the  hole 

o bolt  bearing 

o bolt  shear 

o shear  tear-out 

Tension  at  the  hole  is  frequently  the  dominant  failure  mode.  The  tension 
stress  consists  of  two  components,  a tensile  stress  due  to  the  bolt  bearing 
on  the  hole,  and  a tensile  stress  due  to  a tension  load  running  past  the  hole. 
The  tension  load  running  past  the  hole  is  the  sum  of  the  bolt  loads  upstream. 

I „ Pt  „ pb 

\ °t  = Kti  t(w-d)  + Kt2  te(w-d) 
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= Pt  w Nx 

= observed  stress  concentration  factor  associated 
with  unloaded  open  hole 


= observed  stress  concentration  factor  associated 
with  loaded  hole 

The  bolt  loads  at  a given  row^  and  the  load  running  past  the  row  are 

expressed  as  fractional  parts  and  are  denoted  by  Pb  and  Pt,  respectively. 

In  the  case  of  double-lap  and  single-lap  joints  for  which  an  optimum  design 

4 

can  be  obtained,  a subroutine  based  on  a curve  fit  has  been  coded.  For 
the  stepped- lap  joint,  however,  for  which  only  an  analysis  is  conducted, 
a more  time  consuming  subroutine  employing  matrix  Inversion  is  used.  The 
bolt  load  distribution  is  dependent  on  the  flexibilities  of  the  bolts,  the 
joint  flexibility,  and  the  frictional  force  induced  by  tightening  the  bolts. 
The  frictional  forces,  however,  are  neglected  in  this  investigation  since 
they  depend  on  how  far  the  bolts  were  torqued  at  the  time  of  installation. 

The  observed  stress  concentration  factors  for  composites  are  signifi- 
cantly lower  than  the  elastic,  isotropic  stress  concentration  factors 
which  are  theoretically  derived.  The  observed  stress  concentration  factors 
are  emperically  derived  from  test  data  conducted  for  two  different  composi- 
tions. The  two  compositions  tested  differ  in  the  percentages  of  0°  plies. 


1.  Hart-Smith,  L.  J.,  "Adhesive-Bonded  Double  Joints",  NASA  CR-112235, 
January  1973. 

4.  Yen,  S.  W. , "Investigation  of  Load  Distribution  Among  Fasteners  in  a 
Multiple  Row  Double-Cover  Butt  Joint",  MDC  J5049-01  Douglas  Aircraft 
Company  Report,  1971. 
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37-1/2%  0°-p1ies 

Ktl  = .3923  - .3512  d + (.6077  + .3512  d)  K 
Kt2  = .5949  - .2341  d + (.4051  + .2341  d)  K 


tl 


t2_ 


These  expressions  differ  somewhat  from  those  given  in  reference  1,  due 

to  the  inclusion  of  a size  effect  which  is  introduced  by  the  diameter,  d. 

It  should  be  noted  that  the  above  expressions  are  based  on  very  limited  data, 

especially  with  respect  to  the  size  effect.  As  more  data  becomes  available, 

these  expressions  should  be  updated  accordingly.  The  elastic,  isotropic 

stress  concentration  factors  associated  with  loaded  and  open  holes  are  denoted 

in  the  above  expression  as  and  Kt2  respectively.  They  are  presented 

e e 

in  reference  1 and  are  theoretically  derived. 


Ktl  = 2 + (1  - d/w)' 


t2_ 


= 2 + (w/d  - 1)  - 1.5 

(w/d  + 1) 


6 


8 = 1.5  - .5/(e/w)  for  e/w  < 1 

9 = 1.0  for  e/w  - 1 
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The  effective  thickness  which  is  assumed  to  react  the  bearing  stresses 
at  the  loaded  holes  is  not  allowed  to  exceed  the  bolt  diameter.  Although 


this  demand  is  rather  severe  and  perhaps  overly  conservative,  some  limitation 
is  necessary  to  account  for  bolt  bending  and  a host  of  other  nebulous  effects. 
Mathematically,  the  effective  thickness  is  defined  by 


te  = t for  t-  d 

tg  = d for  t > d 

For  the  case  of  the  unsupported  single-lap  bolted  joint,  there  is  an 
additional  amount  of  tension  at  the  hole,  caused  by  eccentrically  induced 
bending  stresses.  This  bending  stress  is  derived1  for  the  case  of  bonded 
joints  and  is  a maximum  at  the  edges  of  the  bonded  overlap.  The  bending 
stresses1  are  assumed  for  the  case  of  bolted  joints,  but  with  an  overlap 
distance  of  2C  taken  as  the  distance  between  the  outer  rows  of  bolts. 


In  the  case  of  single-lap  joints,  Nx  should  be  limited  to  half  the  value 
of  equation  (5),  which  is  based  on  double  shear.  The  maximum  loading  which 
can  be  applied  before  a shear  tear-out  failure  occurs  at  a particular  bolt 
row  is  determined  from 


F t . 

Nx  = 2(e/d  -.5)  — — — S- 

Pb  W 


. (6) 


where  the  edge  distance  e/d  = 3 is  used.  It  should  be  noted  that  this  formula 
is  valid  only  when  there  are  sufficient  cross  plies  as  in  the  case  for  the 
two  composite  mixes  considered  here.  With  an  insufficient  number  of  cross 

5 

plies,  it  is  pointed  out  that  no  amount  of  edge  distance  will  prevent  shear 
tear-out. 


i Finally,  for  each  bolt  row,  there  are  four  computed  values  of  Nx  computed 

from  equations  (3)  thru  (6),  each  associated  with  a failure  mode.  Denoting 
the  maximum  of  these  four  values  as  N^,  the  margin  of  safety  for  a particular 
bolt  row  is  expressed  in  terms  of  the  applied  joint  load,  Nx,  by 


The  allowable  joint  load  is  then  that  values  of  nx  which  yields  a M.S.  = 0.00 
for  one  of  the  bolt  rows. 


The  joint  weights  which  are  minimized  are  actually  weight  penalties.  The 
weight  penalty  of  a joint  is  the  additional  weight  due  to  the  joint's  splicing 
function.  For  the  single-lap  joint,  it  consists  of  the  weight  associated 
with  one  of  the  thicknesses  extending  over  the  overlap  distance,  less  the  bolt 
holes  thru  both  thicknesses,  plus  the  weight  of  the  bolts  themselves.  For  the 

5.  Hart-Smith,  L.  J.,  "Bolted  Joints  In  Graphite-Epoxy  Composites", 

NASA  CR-1 44899,  1977. 
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stepped  lap  and  bolted  lap  joints,  the  weight  penalty  consists  of  the 
weight  of  the  joint  buildup,  less  the  nominal  thickness,  less  the  bolt 
holes,  plus  the  bolts  themselves.  It  should  be  noted  that  a comparison 
of  a single-lap  joint  weight  with  a double-lap  joint  weight  is  unfair 
since  the  single-lap  joint  is  not  built  up  and  the  advantage  of  reduced 
nominal  stresses  is  not  utilized. 


DESCRIPTION  OF  NUMERICAL  METHODS 

The  search  algorithm  employed  for  determining  an  optimum  joint  design 
is  a direct  brute-force  search  procedure.  The  number  of  bolt  rows  is 
Incremented  from  one  to  four.  The  diameter  is  incremented  in  sixteenths 
from  3/16  to  1-1/2.  For  each  combination  of  bolt  diameter  and  number  of 
bolt  rows,  the  thickness,  t,  and  the  w/d  ratio  are  determined  for  minimum 
joint  weight. 


From  equation  (5) 


N 


x 


Fsu  (bolt)  d 


d 

w 


so  that  an  upper  bound  for  w/d  is 


w/d 


(bolt)  d 
^max 


The  lower  bound  for  w/d  depends  on  whether  a bearing  design  Is  demanded. 
When  that  is  the  case,  from  equations  (3)  and  (4) 

> ^ 

+ 3K  (1  - d/w) 


tFtu  (1  - d/w) 


Pt  Ktl  + Pb  Wtc  Kt2 


d/w  1 


1 

! 


j 
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The  lower  bound  on  w/d  can  be  determined  by  solving  these  two  Inequalities 
for  w/d  and  t simultaneously  with  an  Iterative  process.  That  this  process 
does  Indeed  lead  to  a lower  bound  on  w/d  can  at  least  be  seen  Immediately 
for  other  than  the  unsupported  single-lap  joint.  The  first  Inequality  leads 
to 

"/d"  (\,  Ktl  + t'tc  Kt2> 

In  addition  to  the  above  constraints  Imposed  on  w/d,  the  following 
practical  limits  on  the  feasible  range  of  w/d  have  been  observed. 

3 - w/d  - 12 

With  the  range  on  w/d  thus  defined,  the  value  of  w/d  which  produces  the 
least  weight  joint  is  computed  using  a quadratic  interpolation  routine. 

Thus  for  any  combination  of  bolt  diameter,  d,  and  number  of  bolt  rows,  m, 
the  value  of  w/d  and  t are  determined.  This  then  represents  one  feasible 
design,  and  there  are  as  many  as  88  of  these  feasible  designs.  The  optimum 
design  is  then  taken  as  the  one  with  the  least  joint  weight. 

LIMITATIONS: 

Although  most  of  the  limitations  associated  with  the  math  model  have  already 
been  noted,  they  are  reported  here  for  sake  of  completeness. 
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Double-Lap  Joint 


o maximum  of  4 bolt  rows 

0 3/16  - bolt  dia.  - 1-1/2  (in  1/16  increments) 

o two  compositions  of  graphite  epoxy  material 

25%  0°  plies 
37-1/2%  0°  plies 
o two  bolt  materials 

steel 
titanium 

o same  bolt  sizes  for  all  rows 
o same  bolt  spacing  for  all  rows 

o edge  distance  of  3D 

o row  pitch  of  6D 

o balanced  joint 

o slope  for  joint  buildup  of  15:1 

o 3 - w/d  - 12 


Stepped-Lap  Joint 

o maximum  of  9 bolt  rows 

o 3/16  - bolt  dia.  - 1-1/2  (in  1/16  increments) 
o two  compositions  of  graphite-epoxy  material 
25%  0°  plies 
37-1/2%  0°  plies 


o 


two  bolt  materials 
steel 
titanium 


i 


o bolt  rows  at  center  of  steps 
o slope  for  buildup  of  15:1 
o 3 - w/d  - 12 
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Single-Lap  Joint 


o maximum  of  4 bolt  rows 

o 3/16  - bolt  dla.  - 1-1/2  (In  1/16  Increments) 
o two  compositions  of  graphite-epoxy  material 

25%  0°  plies 
37-1/2%  0°  plies 
o two  bolt  materials 

steel 
titanium 

0 same  bolt  sizes  for  all  rows 
o same  bolt  spacing  for  all  rows 

o edge  distance  of  3D 

o row  pitch  of  6D 

o equal  thicknesses 

o no  buildup 

o 3 - w/d  - 12 


Many  of  the  above  range  restrictions  are  due  to  very  limited  test  data. 
These  restrictions  can  be  reduced  as  more  complete  test  data  becomes  available. 
With  the  above  limitations,  a celling  on  joint  loads  of  about  nx  = 40,000  lb. 
per  Inch  should  be  observed.  Designs  for  greater  joint  loadings  may  become 
questionable  and  impractical. 


Although  many  test  cases  of  bolted  joints  were  run  for  checkout  purposes, 
future  runs  may  indicate  limitations  not  presently  appreciated.  Additional 
limitations  may  conceivably  be  required  for  the  iteration  and  Interpolation 
algorithms  which  are  used,  even  though  these  codes  have  proved  successful  for 
all  the  checkout  cases  run  to  date. 
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The  joint  computer  programs  have  been  run  with  both  IBM  and  CDC  single 
precision  accuracy  with  no  problems. 

SOLUTION  ACCURACY 

The  most  vulnerable  area  to  inaccuracy  of  the  solution  is  the  test 
data  which  formed  the  basis  of  the  observed  stress  concentration  factors. 
For  example,  the  size  effect  in  these  expressions  was  determined  from 
results  of  only  two  bolt  diameters.  Also,  the  weight  of  the  bolts  is 
subject  to  some  error,  depending  on  the  type  of  bolt  and  nut  combination 
used.  The  effective  thickness  used  to  react  bolt  bearing  may  be  a source 
of  considerable  conservation.  The  overall  accuracy  of  the  analytic  solu- 
tions presented  here  are  assumed  to  be  adequate  for  the  four  types  of 
bolted  joints  covered. 
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